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A simple,  approximate,  asymptotic  high  frequency  result  is  developed 
in  this  paper  for  the  electromagnetic  surface  fields  excited  by  infinitesimal 
electric  or  magnetic  current  sources  which  are  placed  on  a smooth,  perfectly- 
conducting  convex  surface  of  any  shape.  These  surface  fields  propagate 
along  Keller's  surface  ray  paths,  and  their  asymptotic  description  is 
uniformly  valid  along  the  ray  including  the  immediate  vicinity  of  the 
source.  Furthermore,  as  the  curvature  of  the  surface  becomes  vanishingly 
small,  this  r 
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a planar,  perfectly-conducting  surface.  The  development  of  this  result 
is  heurlstically  based  upon  the  form  of  the  asymptotic  solutions  for  the 
surface  fields  which  are  first  obtained  for  the  canonical  circular 
cylinder  and  spherical  geometries.  These  solutions  are  then  generalized 
to  the  arbitrary  convex  surface  via  the  local  properties  of  wave 
propagation  at  high  frequencies.  The  effect  of  torsion  associated 
with  the  surface  rays  Is  clearly  Identified  in  this  result  through  the 
presence  of  a factor  T/io  where  T Is  the  surface  ray  torsion  and 
k is  the  surface  curvature  along  the  ray  direction.  Some  currently 
available  numerical  results  based  on  this  analysis  are  indicated 
for  the  circular  cylinder  case.  Included  in  these  numerical  results 
are  mutual  admittance  calculations  for  a pair  of  identical  (axial 
or  circumferential)  slots  in  a circular  cylinder;  these  numerical 
results  show  very  good  agreement  with  the  results  based  on  the  exact 
(modal)  series  solution.  Additional  numerical  calculations  and  experi- 
mental verifications  for  more  general  surface  shapes  are  being  planned 
for  the  near  future. 
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f I.  INTRODUCTION 

I 

I An  approximate  asymptotic  high  frequency  result  is  developed  in  this 

I paper  for  describing  the  electromagnetic  field  in  the  immediate  vicinity 

! of  a smooth,  perfectly-conducting  convex  surface  which  is  excited  either  by 

f an  infinitesimal  electric  current  source  located  in  the  Immediate  vicinity 

of  the  surface,  or  by  an  infinitesimal  magnetic  current  source  located 
directly  on  the  surface.  The  orientation  of  the  electric  current  source  is 
assumed  to  be  along  the  surface  normal;  whereas,  the  orientation  of  the  mag- 
netic current  source  is  assumed  to  be  along  the  surface  tangent.  Thus,  this 
result  essentially  constitutes  an  asymptotic  approximation  to  the  electric 
or  magnetic  type  dyadic  surface  Green's  function  or  its  appropriate  spatial 
derivatives,  for  the  above  orientations  of  the  infinitesimal  electric  and 
magnetic  current  sources.  One  may  therefore  calculate  the  electromagnetic 
field  in  the  immediate  vicinity  of  a smooth,  convex,  conducting  surface 
which  is  excited  either  by  an  aperture,  or  a thin  short  stub  (monopole) 
antenna  located  on  that  surface  by  integrating  this  asymptotic  dyadic  surface 
Green's  function  (or  it's  appropriate  derivatives)  over  the  source  distri- 
bution either  in  the  aperture,  or  along  the  length  of  the  stub,  respectively. 
The  source  distribution  in  these  cases  may  be  defined  in  terms  of  an  equivalent 
tangential  magnetic  current  distribution  in  the  aperture  [1],  or  in  terms 
of  an  equivalent  electric  current  distribution  along  the  stub,  both  of  which 
are  assumed  known.  This  result  for  the  electromagnetic  field  generated  by 
an  aperture  or  stub  antenna  on  a convex,  conducting  surface  may  be  employed 
to  calculate  the  mutual  coupling  between  a pair  of  aperture  or  stub  antennas 
located  on  that  surface  via  a procedure  based  on  the  reaction  theorem  [2]. 

The  details  of  such  an  analysis  of  the  mutual  coupling  between  a pair  of 
radiating  elements  on  a smooth,  convex  conducting  surface  are  presented  in 
this  paper  for  the  case  in  which  one  of  the  radiating  elements  is  an  aperture 
with  the  other  element  being  either  an  aperture,  or  a stub;  and  also  for 
the  case  In  which  both  of  the  radiating  elements  are  stubs.  These  radi- 
ating elements  are  assumed  to  be  short  and  thin  for  reasons  that  will  be 
Indicated  later  In  Section  V;  In  particular,  the  aperture  is  assumed  to  be 
a thin  rectangular  slot  which  is  less  than  a half  wavelength  long,  and  the 
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stub  (or  monopole)  is  assumed  to  be  a thin  wire  which  is  less  than  a quarter 
wavelength  long.  The  problem  of  mutual  coupling  is  of  interest  in  the 
following  way.  A knowledge  of  the  mutual  coupling  between  antennas  on  a 
smooth  convex  surface  is  essential,  for  example,  in  the  design  of  conformal 
arrays,  or  in  estimating  the  electromagnetic  compatibility  of  various 
multi -function  antennas  that  are  located  on  such  a surface.  In  the  present 
context,  the  smooth  convex  surface  could  locally  model  the  region  of  an 
aircraft  or  spacecraft  fuselage  on  which  the  conformal  array,  or  a pair  of 
multi -function  antennas  may  be  located. 

In  the  present  asymptotic  solution  the  surface  field  due  to  an 
infinitesimal  current  source  which  is  placed  on  a convex  surface  is  associated 
with  Keller's  surface  rays  [3,4,5]  which  traverse  geodesic  paths  on  the 
surface  between  the  source  and  the  field  points,  respectively.  In 
particular,  the  field  of  these  surface  rays  is  expressed  in  terms  of  Fock 
type  integrals,  and  this  field  representation  is  uniform  in  the  sense  that 
it  remains  valid  in  the  neighborhood  of  the  source  and  within  the  shadow 
boundary  transition  (or  penumbra)  region;  whereas,  sufficiently  far  from 
the  source  it  can  be  transformed  via  Cauchy's  residue  theorem  into  a rapidly 
converging  series  representation  which  is  commonly  referred  to  as  the 
creeping  wave  modal  series.  Furthermore,  as  the  curvature  of  the  surface 
becomes  vanishingly  small,  this  result  in  terms  of  the  Fock  type  integrals 
recovers  the  known,  exact  solution  for  the  planar  surface.  Basically,  this 
asymptotic  result  is  valid  provided  that  the  surface  curvature  is  neither 
large,  nor  rapidly  varying.  When  the  source  and  the  field  points  are  not 
directly  on  the  surface,  but  are  in  the  immediate  vicinity  of  the  surface, 
the  field  is  expressed  in  terms  of  a Taylor  expansion  which  involves  the 
"soft"  and  "hard"  type*  Fock  integrals  for  the  surface  field  and  its 


♦the  "soft"  and  "hard*1  terminology  employed  here  refers  to  the  soft  (or 
Dirichlet)  and  hard  (or  Neumann)  boundary  conditions  encountered  in 
acoustics.  The  application  of  these  boundary  conditions  leads  to  the 
acoustic  (or  scalar)  wave  solutions  of  the  soft  and  hard  types. 


derivatives.  The  soft  and  hard  type  Fock  integrals  denoted  by  11(5)  and 
vU)  which  appear  in  all  of  the  field  expressions  developed  here  are 
tabulated  by  Logan  [6],  and  their  derivatives  can  be  calculated  numerically 
in  an  accurate  manner.  Thus,  the  present  result  not  only  imparts  a simple 
physical  interpretation  in  terms  of  surface  rays,  but  it  is  also  in  a 
form  which  is  convenient  for  engineering  applications. 


It  is  noted  that  several  other  asymptotic  solutions  for  the  surface 
currents  excited  by  a tangential  magnetic  current  moment  on  a convex 
conducting  surface  have  been  obtained  previously  by  various  authors. 
Hasserjian  and  Ishlmaru  presented  an  asymptotic  solution  pertaining  to  the 
convex  cylinder  geometry  [7];  they  obtained  terms  to  lowest  order  in 

and  where  s is  the  arc  length  of  the  geodesic  ray  path  on  the  surface 

from  the  source  to  the  observation  point,  pg  is  the  radius  of  curvature  of 
the  surface  along  the  ray  direction,  and  k is  the  free  space  wavenumber. 

The  nature  of  their  solution  is  quite  similar  to  that  given  earlier  by 
Wait  [8]  for  a large  conducting  sphere  excited  by  a stub  (or  a short 
monopole)  antenna.  Later,  Hwang  and  Kouyoumjian  [9]  constructed  an 
asymptotic  solution  for  the  convex  cylinder  case  within  the  systematic  frame- 
work of  the  geometrical  theory  of  diffraction  (GTD)  [3,4,5].  In  order  to 
perform  this  extension  of  the  GTD,  they  introduced  the  modal  launching  and 
attachment  coefficients  for  torsional  surface  rays.  The  launching  and 
attachment  coefficients  are  proportional  to  Keller's  surface  ray  modal  dif- 
fraction coefficients  [4,5];  these  launching  and  attachment  coefficients  were 
introduced  earlier  in  [1]  for  the  case  of  non-torsional  surface  rays.  The 
result  developed  in  [9]  has  been  used  in  [10]  for  analyzing  a simple, 
conformal,  cylindrical  array  problem.  In  [9],  the  vector  nature  of  the 
surface  current  (or  the  tangential  surface  magnetic  field)  is  expressed 
rather  compactly  in  a dyadic  form  by  employing  the  unit  tangent  and 
binormal  vectors  fixed  in  the  local  surface  ray  coordinates  at  the  source 
and  the  field  points,  respectively.  As  a result,  this  GTD  formulation  in 
[9]  provides  a simple  physical  interpretation  for  the  manner  in  which  the 
surface  magnetic  field  is  excited  by  a tangential  magnetic  current  source 
on  a convex  cylindrical  surface.  As  the  field  point  approaches  the 


neighborhood  of  the  source,  the  creeping  wave  modal  series  representation;  I 

i.e. , the  GTD  representation  in  terms  of  the  modal  launching  and  attachment 

coefficients  becomes  slowly  convergent,  and  it  is  then  transformed  in  [9] 

into  the  same  Fock  type  integrals  as  those  employed  in  [7,8]  and  in  this 

paper.  Such  a procedure  was  mentioned  earlier  where  it  was  noted  that  these 

Fock  type  integrals  can  be  transformed  via  the  Cauchy  residue  theorem  into 

the  creeping  wave  residue  (or  modal)  series  which  are  rapidly  convergent 

for  large  ks.  In  conformal  array  applications,  the  radiating  elements  are 

not  necessarily  far  apart  so  that  in  the  calculation  of  the  mutual  coupling  !' 

between  two  closely  spaced  array  elements  it  becomes  necessary  to  evaluate 

the  surface  field  for  small  values  of  ks.  In  the  latter  case,  higher  order 

terms  in  1/ks  become  significant  in  the  surface  field  calculation,  and  they  - 

therefore  must  be  included  for  accuracy;  these  higher  order  terms  in  1/ks  are 

not  present  in  [7],  [8]  and  [9].  Also,  the  result  in  [9]  and  [10]  based  on 

the  GTD  indicates  that  the  component  of  the  surface  magnetic  field  along  the 

tangent  to  the  ray  path  may  be  expressed  only  in  terms  of  the  acoustic  soft 

type  surface  ray  field,  and  likewise  the  component  along  the  binormal 

direction  to  the  ray  may  be  expressed  only  in  terms  of  the  acoustic  hard 

type  surface  ray  field;  however,  such  a simplified  field  decomposition 

does  not  appear  to  provide  a sufficiently  complete  field  description  as 

indicated  in  some  of  the  more  recent  analyses  [11,12].  The  GTD  analysis  of  i 

[9]  which  represents  an  asymptotic  approximation  only  to  first  order  ! 

nevertheless  presents  very  useful  and  important  concepts.  In  [11],  Chan  et.  \ 

al.  obtain  asymptotic  expressions  for  the  surface  fields  on  a convex 

cylinder,  and  also  on  a semi-infinite  cone.  The  solution*  in  [11]  is 

expressed  in  the  local  surface  ray  coordinates  as  in  [9],  and  it  contains 

terms  (involving  the  Fock  type  integrals)  to  lowest  order  in  1 /kpg  and  to 

two  orders  in  1 / ( ks ) . As  the  curvature  of  the  convex  surface  becomes 

vanishingly  small,  the  solution  in  [11]  thus  recovers  the  exact  solution 


♦Three  different  solutions  are  presented  in  [11];  here  we  refer  to  their 
solution  which  is  designated  as  Full  Formula  (and  it's  variant  for  the 
circumferential  current  element). 
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for  the  planar  surface  only  to  terms  in  l/(ks)  , but  not  to  terms  in  1 / ( ks ) . 
The  solutions  in  [7],  [9]  and  [11]  for  the  convex  cylindrical  surface  are 
obtained  via  a generalization  of  their  respective  asymptotic  solutions 
which  were  developed  first  for  the  simpler,  canonical  circular  cylinder 
geometry.  Such  a generalization  of  the  canonical  solutions  is  based  on  the 
principle  of  locality  of  wave  propagation  at  high  frequencies  as  employed  in 
the  6TD  procedure.  Further  generalizations  to  conical  surfaces  as  in  [11], 
or  to  arbitrary  convex  surfaces  may  be  performed  heuristically  in  a similar 
fashion.  More  recently,  Lee  [12]  has  presented  an  approximate  asymptotic 
solution  for  an  arbitrary  convex  surface  such  that  it  recovers  the  exact 

O 

planar  result,  i.e. , up  to  terms  in  1/ (ks ) , as  the  surface  curvature  becomes 
vanishingly  small.  The  solution  in  [12]  is  also  expressed  in  terms  of  the 
unit  vectors  fixed  in  the  surface  ray  coordinates  of  [9].  The  solution  in 
[12]  is  obtained  by  first  modifying  an  asymptotic  solution  pertaining  to 
the  canonical  sphere  geometry;  this  "modified"  sphere  solution  is  then  sub- 
sequently generalized  heuristically  to  treat  an  arbitrary  convex  surface 
via  the  principle  of  locality  of  high  frequency  wave  propagation.  In 
particular,  an  additive  type  term  is  introduced  into  the  solution  for  the 
sphere  problem  in  [12]  in  an  ad  hoc  fashion,  to  construct  that  "modified" 
sphere  solution.  As  indicated  in  [12],  the  only  justification  for  including 
this  "ad  hoc"  term  is  that  it  is  essential  for  the  purposes  of  obtaining 
accurate  results  when  the  "modified"  sphere  solution  is  applied  to  the 
circular  cylinder  geometry;  furthermore,  when  this  modified  solution 
is  applied  back  to  the  sphere  geometry,  it  does  not  yield  numerical  values 
which  are  significantly  different  from  those  based  on  the  un-modified 
sphere  solution.  Due  to  the  fact  that  the  solution  for  the  arbitrary 
convex  surface,  and  hence  for  the  convex  cylinder  case  in  [12]  is  based 
on  the  sphere  solution  to  begin  with,  it  therefore  always  contains  both, 
the  acoustic  "hard"  and  "soft"  type  Fock  integrals  corresponding  to  the 
two  types  of  scalar  potentials  from  which  the  sphere  solution  is  con- 
structed. However,  it  is  well  known  that  for  an  axially  directed  magnetic 
current  source  on  a perfectly-conducting  cylinder,  only  the  scalar  potential 
corresponding  to  the  acoustic  "hard"  function  is  required  to  completely 
describe  the  electromagnetic  field  which  is  generated  by  that  source. 
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It  is  noted  that  the  field  solutions  in  [7],  [9], and  C ll3 * contain  only  the 
acoustic  hard  type  potential  for  the  case  of  an  axial  magnetic  current 
source  on  a convex  cylinder.  An  interesting  difference  between  the 
cylinder  result  in  [11]  and  the  result  in  [12]  is  that  a "cross  term" 
exists  in  the  dyadic  representation  for  the  surface  field  of  [11];  whereas, 
it  is  absent  in  [12]  since  the  latter  solution  is  based  on  the  sphere 

O 

solution.  Such  a cross  term  which  involves  higher  order  range  (or  (ks)“  ) 
dependent  terms  is  also  absent  in  [7]  and  [9],  because  the  cylinder 
solutions  in  [7]  and  [9]  do  not  contain  these  higher  order  range  terms. 


In  this  paper,  an  asymptotic  result  for  the  surface  fields  on  an 
arbitrary,  convex,  perfectly-conducting  convex  surface  is  obtained  by 
generalizing  asymptotic  solutions  which  are  first  constructed  in  Appendices 
I and  II  for  the  surface  fields  on  canonical,  perfectly-conducting, 
circular  cylinder  and  spherical  geometries,  respectively.  The  solution 
for  the  sphere  problem  in  Appendix  II  is  essentially  based  on  the  method 
employed  by  Fock  [13].  The  excitation  for  these  canonical  circular  cylinder 
and  sphere  problems  analyzed  in  Appendices  I and  II  includes  not  only  a 
magnetic  current  source,  but  also  an  electric  current  source.  On  the  other 
hand,  only  the  case  corresponding  to  the  magnetic  current  source  excitation 
is  considered  in  [7,9,11,12],  The  canonical  solutions  developed  in  the 
Appendices  are  summarized  in  a dyadic  form  similar  to  that  in  [9,11,12] 
in  terms  of  the  surface  ray  tangent  and  bi normal  vectors.  It  is  noted  that 
the  dyadic  surface  field  representation  for  the  canonical  cylinder  geometry 
contains  a cross  term  similar  to  that  present  in  [11],  but  the  form  of  this 
cross  term  is  not  identical  to  that  in  [11],  Furthermore,  the  solution 
for  the  axial’y  directed  magnetic  current  on  the  cylinder  contains  only 
the  acoustic  hard  type  potential  as  in  [7,9,11].  The  case  of  magnetic 
current  source  excitation  is  summarized  in  Section  II,  and  likewise  the 
electric  current  source  excitation  case  is  summarized  in  Section  III. 

Based  on  the  information  present  in  the  results  for  both,  the  canonical 
cylinder  geometry  and  the  canonical  sphere  geometry,  these  results  are  in 
turn  heuristically  generalized  with  the  aid  of  the  local  properties  of 
high  frequency  wave  propagation,  to  treat  an  arbitrary  convex  surface 


1 
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j 


6 


geometry  as  Indicated  In  Section  IV.  Only  terms  to  lowest  order  In  ^ 

are  retained  In  this  analysis.  The  effect  of  torsion  associated  with 
the  surface  rays  is  clearly  identified  in  this  solution  through  the 
presence  of  a factor  T/k,  where  T is  the  surface  ray  torsion  and 
K*(Pg)"^.  As  the  surface  curvature  becomes  vanishingly  small,  the 
present  solution  recovers  all  terms  in  the  known,  exact  solution 
for  the  planar  geometry.  Expressions  for  mutual  coupling  between  slots 
and/or  stubs  (monopoles)  are  presented  in  Section  V in  terms  of  the 
results  obtained  earlier  in  Section  IV.  Some  preliminary  numerical  cal- 
culations indicating  the  accuracy  of  the  results  developed  in  this 
paper  are  also  shown  in  Section  V. 

II.  SUMMARY  OF  ASYMPTOTIC  SOLUTIONS  TO  THE  CANONICAL 

PROBLEMS  FOR  THE  MAGNETIC  CURRENT  SOURCE  CASE 

As  mentioned  earlier,  the  problems  involving  the  electric  and  magnetic 
current  source  excitation  of  the  circular  cylinder  and  sphere  geometries 
which  have  been  analyzed  in  Appendices  I and  II,  respectively,  serve  as 
canonical  problems  in  the  sense  that  their  solutions  are  instrumental  to 
the  development  of  a solution  for  an  arbitrary  convex  surface,  as  indicated  in 
Section  IV.  Hence,  the  results  obtained  in  Appendices  I and  II  are 
summarized  in  this  section  for  the  case  of  magnetic  current  source  exci- 
tation; whereas,  corresponding  results  for  the  electric  current  source 
excitation  are  summarized  in  Section  III. 

The  results  for  the  circular  cylinder  are  summarized  in  part  A of  this 
section,  and  the  results  for  the  sphere  case  are  summarized  in  part  B. 

These  results  are  expressed  below  in  an  invariant  form  in  terms  of  the 
local  surface  ray  coordinates  as  in  [9]. 

A.  Asymptotic  Results  for  the  Canonical 
Circular  Cylinder  Problem 

The  circular  cylinder  geometry  of  interest  is  illustrated  in  Figure  1. 
The  infinitesimal  magnetic  current  source,  ft  which  excites  the  circular 
cylinder  may  be  represented  in  terms  of  the  Dirac  delta  function  by 
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R-  ■ Pm  5 ( |r-r 1 1). 


(1) 


The  source  R is  assumed  to  be  placed  at  on  the  surface  of  the 
cylinder  and  it  is  oriented  tangential  to  the  surface.  The  quantity  pm 
refers  to  the  strength  of  the  infinitesimal  magnetic  current  (or  current 
moment),  and  F is  an  arbitrary  position  vector;  whereas,  F'  is  the  position 
vector  at  the  source  location  P^  on  the  circular  cylinder. 

First,  expressions  for  the  electric  and  magnetic  fields  are  given 
for  the  case  when  the  field  point  is  located  at  a point  PN  which  is 
on  the  surface.  Next,  an  expression  is  given  for  the  normal  component  of 
the  electric  field  when  the  field  point  is  moved  from  P^  to  the  point  P 
which  lies  above  but  in  the  close  vicinity  of  the  surface. 

From  the  details  of  the  analysis  in  Appendix  I,  it  is  observed  that 
the  field  may  be  interpreted  to  arrive  at  P^  from  the  source  point  at  P^ 
via  a geodesic  ray  propagation  path  on  the  surface.  It  is  noted  that  in 
the  case  of  the  cylinder,  the  geodesic  surface  ray  paths  are  helices.  The 
vector  nature  of  this  field  at  P^  may  be  expressed  compactly  in  terms  of  a 
triad  of  unit  vectors  fixed  in  the  surface  ray  coordinates  at  P^  and  P^. 
These  unit  vectors  are  defined  as  follows.  Let  n and  n'  denote  the  unit 
outward  normal  vectors  to  the  surface  at  PN  and  P^,  respectively;  likewise, 
let  t at  P^  and  t'  at  P^  denote  the  unit  vectors  which  are  oriented  along 
the  forward  tangent  to  the  geodesic  surface  ray  path  from  P^  to  P^.  Then 
the  unit  binomial  vectors  at  P^  and  P^  associated  with  the  geodesic  ray 
path  between  these  points  are  defined  by  b = t x n and  S'  = t'  x n', 
respectively.  These  unit  tangent,  normal,  and  binormal  vectors  constitute 
the  triad  of  unit  vectors  associated  with  a surface  ray.  The  surface  ray 
geodesic  path  between  P^  and  P^,  and  the  associated  triad  of  unit  vectors 
at  these  end  points  are  illustrated  in  Figure  1.  In  addition,  the  pair  of 
unit  vectors  (t-j.T^)  and  (t,,^)  along  the  principal  directions  on  the 
cylinder  surface  at  P^  and  PN  are  also  illustrated  in  Figure  1.  It  is  noted 
that  the  following  relationships  exist  between  the  various  unit  vectors. 
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T 


cos6  = t'  * t -j  = T2*b'  = t * t i = x2*b  , 

sinfi  = t'  *t2  = - (t^  *b' ) = - (t-j-E)  = t*x2. 

Referring  to  (2a;2b),  and  (A-34)  of  Appendix  I,  one  may  now  express 
tangential  magnetic  field  at  P^  on  the  surface  due  to  the  source 
Pj!l  in  terms  of  the  unit  vectors  defined  above  as  follows. 


1 


(2a) 

(2b) 


the 
FT  at 


"tCN^Pn,' 


L_ 

*2<£,ef  fe(1" 


v(0+(f)  is  CuCc)-v(c)3>  + 


+ t't< 


fe1-  fe]  k [“(5)-v(U]l) 


G(ks)  (3) 


where  a = radius  of  the  circular  cylinder, 


k Y -jks 

G(ks)  = “TT  » 

kPn  V3  ? 

5 = — ; with  m = (-*■*■)  and  pn  = a/sin  5. 
pg  d 9 


(4) 

(5) 


The  factor  f is  defined  as  l/sin^fi  in  (A-27b)  of  Appendix  I. 

The  quantity  Yq  is  the  free  space  admittance.  The  quantities  k and  s refer 
to  the  free  space  wavenumber  and  the  length  of  the  surface  ray  geodesic 
path  from  P^  to  PN;  whereas,  pg  denotes  the  radius  of  curvature  of  the 
surface  along  the  ray  (or  t)  direction,  as  indicated  previously  in 
Section  I.  The  result  in  Equation  (3)  represents  along  a given  geodesic 
ray  path  from  P^  to  P^.  In  the  circular  cylinder  problem,  there  are  two, 
dominant,  helical  geodesic  paths  from  P^  to  PN,  and  only  one  such  path  is 
indicated  in  Figure  1;  the  other  geodesic  path  (not  Indicated)  encircles 
the  cylinder  in  the  opposite  sense  with  respect  to  the  one  shown  in  Figure  1. 
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The  total  field  is  then  the  sum  of  the  fields  of  each  of  these  rays.  The 
result  in  Equation  (3)  also  applies  to  the  other  geodesic  path  provided 

A A AAA 

that  the  unit  vectors  t',b',t,b,n,  and  the  quantities  s,  and  c are  now 
associated  with  this  other  ray  path.  In  addition,  infinitely  many 
geodesic  surface  ray  paths  corresponding  to  multiple  encirclements  of  the 
cylinder  (before  reaching  the  field  point)  do  exist;  however,  their  con- 
tribution is  negligible  for  large  ka.  The  functions  u(c)  and  v(c)  refer  to 
the  "soft"  and  "hard"  Fock  type  integrals  which  are  defined  in  Equations 
(A-22a)  and  (A-29)  of  Appendix  I.  As  mentioned  earlier,  Logan  has  provided 
tabulated  values  of  the  Fock  integrals  u(?)  and  v ( c) ; in  addition,  he  has 
also  given  rapidly  convergent  series  expansions  to  accurately  represent  these 
functions  [6],  The  series  expansions  for  uU)  and  v (c)  are  also  indicated 
in  [11,12].  Although  the  series  expansions  in  u(s)  and  v(e)  are  available 
in  [6,11,12],  they  are  presented  in  Appendix  V for  the  sake  of  completeness; 
it  is  noted  that  separate  expansions  exist  for  large  and  small  values  of 
S,  and  these  separate  expansions  provide  numerical  values  which  do  overlap 
for  moderate  values  of  £ as  shown  in  [11].  The  result  in  Equation  (3)  is 
therefore  in  terms  of  functions  which  can  be  computed  rather  easily.  One 
may  express  this  result  in  Equation  (3)  more  compactly  as 

fe-]v(cj)G(ks) 

, (6) 

since  x j t-J  + x£x2  = I$,  where  l$  denotes  an  identity  dyad  on  a surface  so 

that  p •!  sp^,  and  (t£*6')6  + (x2*t‘)t  = x2.  Even  though  the  result  in 
Equation  (6)  is  more  compact  than  the  one  in  Equation  (3),  the  expression 
in  Equation  (3)  indicates  more  directly  the  fact  that  the  field  fl!j!  which  is 
produced  by  the  axial  or  the  x^-component  of  the  source  pm  is  substantially 
different  from  that  which  is  produced  by  the  circumferential  or  the  x£- 
component  of  pm.  In  particular,  the  x^  or  the  axial  component  of  p^ 
generates  a field  h!£  which  contains  only  the  hard  type  Fock  Integral  v(0 


«;(PNhpra-(6'b^  fen-  kj  v(s)+t't[ji(i- 

* V;2;2  (f)  fet“(«)-*(«)])s(ks) 


as  expected;  whereas,  the  -ri-component  of  p generates  both,  u(?)  and  I 

“ III  * 

v(0  type  terms.  That  the  component  of  generates  only  a vU)  type  ] 

term  Is  in  agreement  with  the  results  in  [7,9,11];  on  the  other  hand,  the  1 

result  in  [12]  indicates  that  this  axial  component  of  pm  generates  both,  | 

u(0  and  v(5)  type  functions  in  the  cylinder  problem.  Also,  the  "cross"  j 

term  of  the  type  — )G(ks)  in  Equation  (6)  is  also  present  in  [11],  f 

except  that  the  functional  form  of  this  term  in  [11]  is  somewhat  different  j 

from  that  in  Equation  (6).  Such  a "cross  term"  is  not  present  in  the  j 

result  of  [12].  Furthermore,  except  for  the  important  difference  which  j 

exists  in  the  functional  form  of  this  "cross  term",  the  result  given  in  [11]  < 

otherwise  agrees  with  that  of  Equation  (6)  up  to  terms  in  1/ks  and  4 i 

2 ^ t 

l/(ks)  . It  is  noted  that  terms  of  order  1 / (ks )°  are  present  in  Equation  | 

(6),  but  they  are  not  present  in  [11].  | 

j 

From  Equation  (A-36b)  of  Appendix  I,  and  the  relationships  in  Equation  j 

(2a;b),  one  may  similarly  express  the  component  of  the  electric  field  Em  i 

n ? 

which  is  normal  to  the  surface  at  P^,  and  which  is  excited  by  the  source  R" 

at  Pj||,  in  terms  of  the  surface  ray  coordinates  as  follows.  | 

I 

cos«nf  jL[u(0-v(0])s(ks).  (7)  ! 

i 

I 

As  before.  Equation  (7)  represents  the  field  of  a "typical"  surface  ray  path  J 

from  Pj^  to  P^;  there  are  of  course,  two  such  dominant  ray  paths  and  the  field  \ 

of  both  of  these  rays  must  be  included  to  obtain  the  total  field.  j 

j 

The  results  in  Equation  (7)  may  now  be  generalized  such  that  the  field  | 

e!J  is  evaluated  not  on  the  surface  at  PN,  but  is  instead  evaluated  at  a $ 

point  P which  lies  above  the  surface.  It  is  assumed  that  the  point  P lies  | ] 

in  the  immediate  neighborhood  of  the  surface.  In  particular,  the  point  ] 

PN  on  the  surface  represents  the  projection  of  the  point  P along  the  normal  j 

n ■ j 

to  the  surface,  and  the  height  of  the  point  P above  P^  is  denoted  by  d2  ; 

as  shown  in  Figures  1 and  2.  Thus,  i 
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P * p'(  o + d,  ,0,0) 

P 8 P (a  + d2,^, z ) 
r,' «2  AT  Pi;£,*$  AT  P 
f2*^  AT  Pn;  %•$  AT  P, 
b *txn  ; b'«  t'xn' 
t - z « cos 8 - t‘.  z' 


lp,4>) 

‘GEODESIC  PATH 
OF  LENGTH  S 


Figure  1.  Circular  cylinder  geometry. 


\i.p 


Figure  2.  Cross  sectional  view  indicating  the  heights 
d]  and  62  of  the  points  p'  and  P above  the 
surface. 


TTF  = \*tf\n  = d2n. 


(8) 


From  the  comments  at  the  end  of  Appendix  I,  it  is  seen  that  the  results  for 
E^(P)  is  the  same  as  that  in  Equation  (7)  with  the  exception  that  the 
functions  uU)  and  v(0  must  now  be  replaced  by  the  more  general  Fock 
functions  F U,y.|,y2)  and  Fh(c,y-|  ,y2) , respectively,  which  are  defined 
in  Equation  (A-139)  and  Equation  (A-138)  of  Appendix  V.  Hence, 

l>MYo)-Vm{6'S|T-i^Fh({>0,y2>+;$  COSS  Sf  jL{FsU,0,y2)-Fh({.0,y2^. 

•G(ks) 


The  quantities  y^  and  y2  can  be  shown  to  be  related  to  d^  and  d2  as  follows. 

y1  = m”1kd1  = 0 (since  M is  at  P^);  y2  = m-1kd2  . (9b;9c) 

Since  kd2  is  assumed  to  be  much  smaller  than  ka  in  the  present  case  (i.e., 

P is  in  the  immediate  neighborhood  of  the  surface),  a Taylor  series  approxi- 
mation for  the  functions  Fs(c,y-|  ,y2)  and  FhU,y^,y2)  may  ^ ^P10^  as 
indicated  in  Equations  (A-146)  and  (A-145)  of  Appendix  V.  Such  an 
approximation  appears  to  be  valid  only  if  the  value  of  K is  not  too  close  to 
zero;  however,  it  has  the  advantage  that  for  sufficiently  small  ky-j  and  ky2, 

the  functions  F$,  and  F^  may  be  evaluated  in  terms  of  the  derivatives  of 


the  functions 


-j  — — it  u 

Ur 


(C),  and  lj|v(c)»  respectively;  these  derivatives 


with  respect  to  the  argument  c may  be  easily  computed  numerically  in  an 
accurate  fashion;  on  the  other  hand,  accurate  approximations  for  these 
functions  and  their  derivatives  are  available  (see  Appendix  V). 


B.  Asymptotic  Results  for  the  Canonical 
Spnere  Problem  ~~ 


The  sphere  geometry  is  illustrated  in  Figure  3.  The  infinitesimal 
magnetic  current  source  FT  is  located  at  P^  on  the  sphere  and  it  has  been 
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GEODESIC  PATH 
OF  LENGTH  ,* 


P‘«  P'(<H-d,  ,0,0) 
P-  P la+d 2,8,+) 
t * q6 


Figure  3.  Sphere  geometry. 


defined  previously  in  Equation  (1).  Following  the  development  in  Part  A 
for  the  circular  cylinder,  the  tangential  magnetic  field  at  PN  on  a 
perfectly-conducting  spherical  surface  which  is  excited  by  FT  may  also  be 
expressed  in  terms  of  the  local , geodesic  surface  ray  coordinates  as 
follows.  Firstly,  it  is  noted  that  the  geodesic  surface  ray  paths  on  the 
sphere  are  great  circles.  The  triad  of  unit  vectors,  (t'.n'.b')  at  P^|, 
and  (t,n,b)  at  PN,  which  are  associated  with  a geodesic  (great  circle) 
path  that  connects  these  points  are  illustrated  in  Figure  3.  Secondly, 
the  following  relationships  exist  between  the  unit  vector  triad  associated 
with  a surface  ray  path,  and  the  usual  triad  of  unit  vectors  (r,e,$) 
associated  with  the  spherical  coordinates. 


e = t ; <}>  = - b 
A1  so, 

Pm  * t'  - cos  * 


r = n . 

pL  • b'  = - sin  $ . 
m 


0<e<Tr 


(10a);(10b);(10c) 
(lOd); (lOe) 


14 


When  Equations  (lOa-lOe)  above  are  Incorporated  Into  Equations  (A-104a;b;c) 
of  Appendix  II,  one  obtains 

fee-  |f)u(o])D-s(ks). 

(ID 

where 


i 


a quantity  which  is  related  to  the  surface  ray 
divergence  factor, 


and 


(12) 


31 


£ = j2-  , with  s = ae.  (13) 

As  in  Equation  (6),  the  above  result  in  Equation  (11)  is  associated  with 
only  one  geodesic  (great  circle)  path  between  P^  and  PN  as  shown  in 
Figure  3.  A second  geodesic  (great  circle)  path  also  exists  between  P^ 
and  P^;  this  path  is  not  shown  in  Figure  3,  and  it  corresponds  to  the  one 
which  arrives  at  P^  after  being  launched  at  P^  in  a direction  opposite  to 
the  path  shown  in  Figure  3.  The  result  in  Equation  (11)  applies  also  to 
, the  other  geodesic  path  provided  the  unit  vectors,  t',b',t,b,  and  n, 

and  the  quantities,  s,  0,  D and  £ are  now  associated  with  this  other  ray 

path.  The  contribution  from  additional,  multiply  encircling  surface  rays 

2 

is  negligible  for  large  ka.  Except  for  the  factor  D which  appears  with 
the  square  brackets  in  Equation  (11),  this  result  agrees  with  the  one  obtained 
by  Lee  [12].  It  appears  that  D has  been  approximated  by  unity  in  [12].  It 
f is  noted  that  e = it,  or  ks  = ka*  is  the  position  of  a focal  caustic  of  the 

surface  rays  on  a sphere.  At  e = it,  D is  unbounded,  and  the  result  in 
Equation  (11)  must  be  modified  as  indicated  in  Appendix  III.  In  particular, 
the  total  surface  magnetic  field,  p^(tota|).  the  sum  0f  the  fields 

associated  with  "both"  of  the  great  circle  geodesic  paths  from  P^  to  P^ 
must  be  replaced  by 


! 
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Pr  in  the  close 
neighborhood 
of  e=7r 


(14) 


where  the  caustic  correction  factors  Tc^  and  T ^ are  given  in  Equations  (A-128) 
and  (A-130)  of  Appendix  III.  Away  from  6=tt,  the  result  in  Equation  (14) 
accurately  approximates  the  sum  of  the  field  associated  with  each  of  the 
two  geodesic  (great  circle)  ray  paths  from  to  P^  in  which  the  field  of 
either  one  of  these  rays  is  typically  given  by  Equation  (11). 


From  Equation  (10),  and  Equation  (A-109)  of  Appendix  II,  one  may  also 
write  the  normal  component  of  the  electric  field  at  PfJ  due  to  the 
source  FT  at  P^,  in  terms  of  the  surface  ray  coordinates  associated  with  a 
typical  geodesic  surface  ray  path  between  P^  and  PN  as  follows. 

felvt<:>)° G(ks)-  (,5) 

Actually,  the  total  field  at  P^  is  the  sum  of  the  fields  associated  with 
the  two  dominant  geodesic  (great  circle)  paths  which  exist  between  the 
points  P^  and  P on  the  sphere;  the  form  of  the  field  along  each  of  these 
paths  is  as  in  Equation  (15). 


The  result  in  Equation  (15)  may  be  generalized  to  include  the  case 
in  which  the  field  point  is  no  longer  on  the  surface  at  P^,  but  is  moved 
to  a point  P which  lies  above  the  surface;  it  is  assumed  in  this  generali- 
zation however,  that  P is  in  the  close  neighborhood  of  the  surface.  As 
before,  PN  is  the  normal  projection  on  the  surface  of  the  Point  P which 
is  located  a distance,  dg  above  the  surface.  The  relationship  between 
PN  and  P ts  as  given  in  Equation  (8).  From  the  discussion  at  the  end  of 
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Appendix  II,  it  is  seen  that  when  d^O,  the  v(£)  in  Equation  (15)  must  be 
replaced  by  Fh(£,y.|  ,y2)  of  Equation  (A-138)  with  y^O.  The  quantities  y-j 
and  y2  are  as  in  Equations  (9b)  and  (9c),  respectively.  Since  kd2  is 
assumed  to  be  much  smaller  than  ka,  may  aPPr0Xlmated  by  a 

Taylor  series  for  small  y2  as  in  Equation  (A-145)  provided  that  £ is  not 
allowed  to  approach  too  close  to  zero.  As  0->tt,  [>►«,  and  the  total  surface 
electric  field  ^(tota^)  . i.e.,  the  sum  of  the  fields  ^ associated  with 
"both"  of  the  great  circle  geodesic  paths  from  P^  to  PN  must  be  replaced  by 


cm(total ) 
n 


0%ir 


(16) 


where  the  caustic  correction  factor  Tc2  is  given  in  Equation  (A-129)  of 
Appendix  III.  If  P^  is  changed  to  P in  Equation  (16),  then  v(£)  must  be 
replaced  by  ^U.y^yg)  with  y^O. 


III.  SUMMARY  OF  ASYMPTOTIC  SOLUTIONS  TO  THE 
CANONICAL  PROBLEMS  FOR  THE  ELECTRIC 
CURRENT  SOURCE  CASE 

In  this  section,  asymptotic  solutions  for  the  fields  in  the  immediate 
vicinity  of  the  canonical  perfectly-conducting  circular  cylindrical,  and 
spherical  geometries  are  summarized  for  the  case  when  they  are  excited  by 
an  infinitesimal  electric  current  source.  The  electric  current  source  is 
placed  in  the  close  neighborhood  of  the  cylindrical  and  spherical  geometries, 
and  the  orientation  of  this  source  is  assumed  to  be  along  the  surface 
normal.  As  mentioned  previously,  the  asymptotic  solutions  to  these 
canonical  problems  are  developed  in  detail  in  Appendices  I and  II  for  the 
circular  cylindrical  and  spherical  geometries,  respectively.  The  results 
for  the  circular  cylinder  case  are  presented  in  Part  A of  this  section, 
and  corresponding  results  for  the  spherical  case  are  presented 
subsequently  in  part  B. 
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iA.  Asymptotic  Results  for  the  Canonical 

Circular  Cylinder  Problem 

\ 

l The  circular  cylinder  geometry  is  illustrated  as  before  in  Figure  1. 

| The  source  and  the  field  points  are  initially  chosen  to  lie  at  and  PN 

I on  the  surface.  The  electric  current  source  T is  defined  as 

! T«  pe  «(|r-r'|).  (17) 

It  is  noted  that  F1  is  the  position  vector  of  the  source  at  P^  on  the 
cylinder  surface,  and  the  quantity  pg  represents  the  strength  of  the 
infinitesimal  electric  current  (or  current  moment).  From  Equation 
(A-39)  of  Appendix  I,  the  tangential  magnetic  field  at  P^  due  to  7 at 
P|(l  propagates  along  a typical  geodesic  surface  ray  path  as  shown  in 
j Figure  1.  This  field  may  be  expressed  in  terms  of  the  unit  vectors 

t',n',b'  and  t,n,b  fixed  in  the  surface  ray  as, 

flJ(PN)'v-(Y0)”1pe*^,b[u  jLjv(c)+n'  cos6  t2  f £g[uU)-v(0^ G(ks). 

(18) 

There  are  of  course  two  dominant  geodesic  surface  ray  paths  as  mentioned 

i 

earlier;  only  one  such  path  is  shown  in  Figure  1.  The  field  of  either  of 
these  rays  is  given  by  Equation  (18). 

If  the  source  at  P/.  is  now  raised  a distance  d.  above  the  surface 
N 1 

f to  the  point  P'  as  shown  in  Figure  1,  where  P'  is  in  the  close  vicinity  of 
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(19)  since  PN  is  on  the  surface.  Only  y1  * nf'kdj  may  be  non-zero.  Further- 
more, PjJjF'  = |P^P'  |n'  = d^'.  As  before,  the  FsU,y1#0)  functions  may  be 

approximated  by  a Taylor  series  valid  for  sufficiently  small  kd1  (kd1«ka) 
and  for  ? not  too  close  to  zero,  as  indicated  in  Equations  (A-145)  and 
(A-146)  of  Appendix  V. 

From  Equation  (A-42)  of  Appendix  I,  the  component  of  the  electric 
field  which  is  normal  to  the  surface  at  P^,  due  to  T at  P^,  may  be 
expressed  in  terms  of  the  surface  ray  coordinates  associated  with  a 
"typical"  geodesic  surface  ray  path  between  P^  and  P^  as  follows. 

^(PN)'-(V0)-2Pe-S'^[-  fco-  y v(5)+  f jL  [u(5)-v(«)> 

- jL-(l-  y[u(e)-v(c)]jG(ks).  (20a) 


^(PNh-(Yo)"2Pe*n'^v(0- ^0- ij)u(5)+  f ^tu(5)-v(?)^G(ks)  (20b) 

One  notes  that  if  source  and  observation  points  at  P^  and  PN  are  both  moved 
to  the  points  P1  and  P,  respectively  which  lie  above  the  surface  such  that 
= d-j n * with  kd1«ka,  and  FjjF  = d2n  with  kd2«ka,  then  the  u(c)  and 
v(0  in  Equation  (20)  must  be  replaced  by  FsU,y.|,y2)  and  ^^*^1*^2^’ 
respectively.  Thus, 

^(P)'-(Vor2Pe-S'^  fe(l-  4s3Fh(e’yl'y2)+  f l$FsU,yl,y2)  ■ 

- Fh(5.y,.yz)]-  k)  [fs ( « ,y,  .y2)-F„( t ,y, ,y2) ]) G (ks ) , 


with  y^  * m”  kd^. 


B.  Asymptotic  Results  for  the  Canonical 
Sphere  Problem 


The  geometry  of  the  sphere  problem  is  illustrated  in  Figure  3.  The 
source  T as  defined  in  Equation  (17)  is  initially  located  at  P^  on  the 
spherical  surface.  From  Equations  (lOa-lOe),  and  (A-110)  of  Appendix  II, 
the  tangential  magnetic  field  at  on  the  sphere  due  to  the  source  T 
at  PjJ|  may  be  expressed  in  terms  of  the  surface  ray  coordinates  associated 
with  a typical  geodesic  great  circle  ray  path  between  P^  and  PN  as 
indicated  below. 


^(^^(ToJ'Ve-Kb  l-  fe]v(c)jD  G(ks)  . 


As  before,  the  total  field  at  PN  is  the  sum  of  the  fields  associated  with 
the  two  dominant  geodesic  surface  ray  paths  which  exist  between  the  points 
P^  and  PN  on  the  sphere;  the  form  of  the  field  along  each  of  these  paths  is 
given  by  Equation  (22).  As  e-^,  PN  approaches  the  focal  caustic  of  the 
surface  rays;  hence,  the  total  surface  magnetic  field  i.e. , the 

sum  of  the  fields  associated  with  "both"  of  the  great  circle  surface  ray 
paths  from  P^  to  P^  must  be  replaced  by: 

e(total  , f A i-  . -»  \ 

"t  <PN>  ^(V0>'V(n'*|l-yv<5))Tc2  • (23> 

0%r 

The  factor  T g is  defined  in  Equation  (A-129)  of  Appendix  III.  Also,  if 
Pj^  is  moved  to  P'  with  P^P 1 =d^ n ' , and  kd^«ka  as  before,  then  v(c)  in 
Equation  (22)  must  be  replaced  by  Fh(£,y.|,y2)  with  y2=0. 

From  Equations  (lOa-lOe),  and  (A-112)  of  Appendix  II,  the  normal 

— b 

component  of  the  electric  field,  En  at  P^  due  to  J at  P^  is  expressed  as 
follows 


•s's(E  fe('-  t3v‘5>)D  G<ks>' 
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Again,  the  result  in  Equation  (24)  represents  the  field  associated  with  a 
"typical"  geodesic  surface  ray  path  between  the  points  P^  and  PN  on  a 
sphere.  If  the  source  and  observation  points  at  P^  and  P^  are  both  moved 
to  the  points  P'  and  P,  respectively  which  lie  above  the  surface  such  that 
PjJjP^d^n*  with  kd^«ka,  and  PfjP^n  with  kdg^ka,  then  v(c)  in  Equation 
(24)  must  be  replaced  by  F^S.y-j^).  Finally,  if  the  re|ult  in 
Equation  (24)  becomes  singular,  and  the  total  electric  field  Fn^tota^;  i.e. , 
the  sum  of  the  fields  associated  with  "both"  of  the  great  circle  paths 
between  P^  and  PN  must  be  replaced  by 

fe  (5-  k|vU)Tci‘  (25) 

where  the  caustic  correction  factors  T^  is  as  given  by  Equation  (A-128) 
of  Appendix  III. 

IV.  GENERALIZATION  TO  ARBITRARY  CONVEX  SURFACES 

In  this  section,  the  results  which  were  presented  in  Sections  II  and 
III  for  the  canonical  circular  cylinder  and  spherical  geometries  will  be 
generalized  to  the  arbitrary  convex  surface  case.  In  particular,  the 
results  pertaining  to  the  surface  fields  which  are  excited  by  infinitesimal 
electric  or  magnetic  current  sources  located  on  a smooth,  perfectly- 
conducting  convex  surface  of  arbitrary  shape  will  be  developed  in  this 
section. 

The  generalization  of  the  canonical  solutions  to  treat  the  arbitrary 
convex  surface  case  is  based  on  the  locality  of  high  frequency  wave  propa- 
gation as  employed  in  the  GTD  procedure.  Thus,  the  results  pertaining  to 
the  canonical  circular  cylinder  geometry  which  were  presented  in  Sections 
II-A  and  III-A  for  the  magnetic  and  electric  source  excitations,  respectively 
may  be  generalized  via  the  GTD  procedure  to  treat  the  arbitrary  convex 


cylinder  with  a slowly  varying  curvature,  by  assuming  that  the  neighborhood 
of  each  Doint  on  the  convex  cylinder  can  be  aDDroximated  bv  a circular 


! 
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cylinder  of  the  same  local  radius  of  curvature  as  the  convex  cylinder. 

While  such  a generalization  to  the  convex  cylindrical  surface  of  variable 
curvature  is  fairly  direct,  an  analogous  generalization  presented  here  to 
treat  a smooth  convex  surface  of  any  shape  is  somewhat  less  direct,  and 
more  heuristic  in  nature  for  the  following  reasons.  The  circular  cylinder 
constitutes  a "developable"  surface;  thus,  it  possesses  an  infinite  radius 
of  curvature  along  its  "generator"  which  also  constitutes  one  of  the 
"principal  directions"  on  the  surface.  In  addition,  the  surface  geodesics 
are  always  straight  lines  on  the  "developed"  surface.  The  arbitrary  convex 
cylinder  and  the  cone  also  belong  to  the  class  of  developable  surfaces; 
hence,  they  possess  all  of  the  properties  of  developable  surfaces  mentioned 
above  in  connection  with  the  circular  cylinder  geometry.  As  a result  of 
these  common  geometrical  properties  which  characterize  the  cylinder  and 
cone,  it  is  therefore  reasonable  to  conjecture  that  a GTD  type  generalization 
of  the  results  pertaining  to  the  canonical  circular  cylinder  geometry  would 
be  accurate  for  treating  not  only  the  convex  cylinder,  but  also  the  cone 
geometry.  This  conjecture  has  of  course  been  verified  by  Keller  [4]  in  j 

the  case  of  the  GTD  analysis  of  the  scattering  of  waves  by  a smooth  convex 
cylinder  wherein  the  GTD  solution  based  on  the  generalization  of  the  : 

corresponding  canonical  circular  cylinder  solution  was  found  to  agree  with 
the  rigorous  asymptotic  high  frequency  solutions  to  the  problem  of  scattering 
of  waves  by  elliptic,  and  parabolic  cylinder  geometries.  Recently,  Chan 
et.  al.  [11]  have  indicated  that  they  are  indeed  able  to  verify  that  the 
GTD  prescription  for  the  surface  ray  field  on  a cone  which  is  based  on  a 
f generalization  of  the  canonical  circular  cylinder  result,  is  in  agreement 

| with  the  corresponding  rigorous  asymptotic  high  frequency  solution  for  the 

j surface  field  on  a cone;  this  verification  in  [11]  pertains  to  a cone 

[ with  an  acoustic  hard  (or  the  Neumann)  boundary  condition.  On  the  other 

hand,  a convex  surface  with  an  arbitrary  shape  is  in  general  neither  a 
surface  of  revolution,  nor  a developable  surface;  hence,  a generalization 
| of  only  the  results  for  the  canonical  circular  cylinder  geometry  is 

expected  to  be  inaccurate  to  treat  an  arbitrary  convex  surface;  clearly, 
additional  information  is  required  in  order  to  construct  an  asymptotic 
result  which  is  useful  for  treating  this  case  since  one  notes  that  the  two 
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principal  radii  of  curvatures  on  an  arbitrary  convex  surface  are  tn  general 
different  and  finite.  In  addition,  the  torsion  associated  with  the  geodesics 
on  an  arbitrary  convex  surface  is  generally  not  a constant.  In  contrast, 
the  canonical  circular  cylinder  geometry  possesses  a constant  principal 
radius  of  curvature  along  the  circumferential  direction  on  the  surface; 
whereas,  it  exhibits  an  infinite  radius  of  curvature  along  the  other  (or 
the  axial)  principal  direction.  Furthermore,  the  torsion  is  constant  along 
a given  geodesic  path  on  a circular  cylinder.  Of  course,  a simple  canonical 
geometry  for  which  both  principal  radii  of  curvatures  are  finite  is  the 
sphere.  The  canonical  problem  of  the  radiation  by  an  infinitesimal 
electric  or  a magnetic  current  source  on  a perfectly-conducting  sphere  has 
been  treated  in  Appendix  II  as  mentioned  earlier.  However,  the  sphere  has 
a constant  curvature,  and  hence  the  surface  geodesics  on  a sphere  are  great 
circle  paths  which  possess  no  torsion.  A third  canonical  problem  for  which 
a well  known,  exact  solution  is  available  is  the  problem  of  the  radiation 
by  an  infinitesimal  electric  or  a magnetic  current  source  on  a perfectly- 
conducting  planar  surface  of  infinite  extent.  The  planar  geometry  has  zero 
curvature;  and  the  surface  geodesics  for  this  case  are  straight  lines  which 
are  obviously  torsionless.  Even  though  the  convex  cylinder,  and  the 
canonical  sphere  geometries  do  not  possess  all  of  the  general  properties 
associated  with  a convex  surface  of  an  arbitrary  shape,  the  asymptotic 
solutions  pertaining  to  these  geometries  including  the  exact  solution  for 
the  planar  case  nevertheless  do  provide  some  useful  information  based  on 
which  an  approximate  asymptotic  solution  may  be  constructed  for  the 
arbitrary  smooth  convex  surface.  Thus,  an  asymptotic  solution  for  the  surface 
fields  excited  by  infinitesimal  electric  or  magnetic  currents  located  on  a 
smooth,  perfectly-conducting  convex  surface  of  any  shape  is  developed  here 
on  the  basis  of  certain  assumptions.  Some  of  these  assumptions  are  com- 
pletely heuristic  in  nature.  In  particular,  the  development  of  this 
asymptotic  solution  is  based  on  the  following  properties  and  assumptions: 


(a)  An  asymptotic  solution  for  the  arbitrary  convex  surface  should 
be  expected  to  recover  the  solutions  pertaining  to  the  canonical  circular 
cylinder  and  sphere  geometries  when  it  is  specialized  to  the  latter  cases, 
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respectively.  This  solution  should  also  be  expected  to  reduce  to  the  well 
known  result  for  the  planar  surface  in  the  limit  as  the  curvature  of  the 
convex  surface  becomes  vanishingly  small. 

(b)  Since  the  sphere  geometry  exhibits  uniform  (or  isotropic)  surface 
properties  along  different  directions  on  the  surface  that  correspond  to 
the  different  geodesic  paths  which  originate  from  a source  on  the  sphere,  the 
scalar  components  of  the  dyadic  transfer  function  relating  the  surface  field 
at  the  receiving  (or  observation)  point  to  the  launching  (or  source)  point 
are  also  independent  of  these  directions  as  seen  from  the  canonical  sphere 
results  of  Sections  II-B  and  III-B,  respectively.  On  the  other  hand,  the 
circular  cylinder  exhibits  non-uniform  (or  non-isotropic)  surface  properties 
along  the  different  directions  that  again  correspond  to  the  different 
geodesic  paths  originating  from  a source  on  this  surface.  As  a result,  the 
scalar  components  of  the  dyadic  transfer  function  for  the  surface  field  on 
a circular  cylinder  In  general  are  a function  of  the  launching/receiving 
direction  on  the  surface.  It  is  easily  seen  from  the  canonical  circular 
cylinder  geometry  that  this  non-uniform  surface  property  results  directly 
from  the  variation  of  the  geodesic  radius  of  curvature,  pg  as  the 
launching/receiving  direction  of  this  geodesic  path  is  changed  at  the 
source/ receiver  location.  Thus,  as  the  launching  direction  of  the  geodesic 
ray  path  is  changed  at  the  source  (launch)  point  on  a circular  cylinder, 
from  the  principal  direction  along  the  circumference  to  the  other  principal 
direction  along  the  generator  (axial  direction),  the  radius  of  curvature 
Pg  In  the  direction  of  the  geodesic  path  increases  monotonically  from 
Pg=a  to  p -*•.  As  a result  of  the  change  in  pg,  the  torsion  associated 
with  the  geodesic  path  also  increases  from  zero  to  a maximum  value  as  pg 
increases  from  a to  2a.  As  expected,  however,  the  torsion  then  decreases 
back  to  zero  as  pg  increases  from  2a  to  infinity,  since  the  geodesics  along 
the  principal  directions  on  a convex  surface  must  be  torsionless.  Referring 
to  the  results  for  the  surface  fields  on  a circular  cylinder  as  given  in 
Sections  II-A  and  III-A,  one  notes  that  the  dependence  of  the  scalar  com- 
ponents of  the  surface  field  dyadic  transfer  function  upon  variations  in 
pg,  resulting  from  a change  in  the  launching  direction  of  the  geodesic 


path  at  the  source  can  be  Identified  more  definitely  If  the  unit  vectors 
and  In  Equations  (6),  (7),  (18),  and  (20a)  are  expressed  In  terms  of  the 
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unit  vectors  t* , b‘,  t and  b fixed  In  the  surface  rays.  The  corresponding 
results  In  Equations  (11),  (15),  (22),  and  (24)  for  the  sphere  case  have  been 
completely  expressed  In  the  surface  ray  coordinates.  Once  the  cylinder  re- 
sults are  also  completely  expressed  In  an  Invariant  form  In  terms  of  the  same 
surface  ray  coordinates,  the  precise  differences  between  the  cylinder  and 
sphere  solutions  should  become  apparent.  Since  the  goedeslc  surface  ray 
paths  In  general  posses  non-zero  torsion  in  the  case  of  the  cylinder;  whereas, 
they  are  torsionless  in  the  case  of  the  sphere,  it  is  therefore  also  reason- 
able to  expect  that  any  differences  between  the  cylinder  and  sphere  solutions 
may  be  primarily  caused  by  the  effects  of  torsion  associated  with  the  surface 
rays  on  a cylinder. 

In  order  to  express  Equations  (6),  (7),  (18),  and  (20a)  completely  in 
terms  of  the  surface  ray  coordinates,  one  makes  use  of  the  following  relation- 
ships. 


xj  = cosst'-sinsb' 

; t = cos6t-sin6b 

(26a ;26b) 

* sinst'+cosfib' 

; t2  = sinfit+cosfib 

(26c;26d) 

t2t2  = s^n26T,t+cos26b,b+sin6cos6(t'b+b't) 

( 26e) 

I Incorporating  Equation  (26e)  into  Equation  (6)  yields 

[ ^<PN>'\'(i'^[<1-t>',(5>+D2(fe,Zv(5>4-Tofe{u(5)-vU)}]  + 

+ t'tjp2  vU)+jjj  uU)-2(jL)  v(5)J  + 

+ [t'b+b't]T0  jL  {u(0-v(5)})d  G(ks)  . (27) 

i 

f The  factor  D which  is  related  to  the  surface  ray  divergence  factor  is  unity 

I for  the  circular  cylinder;  however,  it  is  introduced  symbolically  in  Equation 

(27)  above  to  facilitate  a direct  comparison  with  the  sphere  results.  The 
factor  T (=cot6)  is  identified  as  a ratio  of  the  surface  ray  torsion  T and 
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the  surface  curvature  K(al/pq)  along  the  ray  direction.  It  is  easily  verified 
that  for  a given  helical  geodesic  surface  ray  path  on  a circular  cylinder, 


Tq  = cots  = T/k 


sin26 

2a 


‘ = (pg> 


sin  <5  # ( 28a  ;28b ;28c) 
a 


Thus,  the  quantity  TQ  is  a constant  for  a given  helical  surface  ray  on  the 
cylinder,  and  indeed,  it  serves  to  uniquely  specify  that  helical  geodesic 
path.  When  T (and  hence  T ) is  allowed  to  vanish  as  in  the  case  of  the  sphere, 
the  above  result  in  Equation  (27)  for  the  cylinder  reduces  exactly  to  the 
result  in  Equation  (11)  for  the  sphere,  except  for  the  ()jj)^v(s)G(ks)  term 
which  is  different  from  the  (j|y)^u(c)G(ks)  present  in  the  sphere  case.  As 
expected,  the  significant  differences  between  Equations  (27)  and  (11)  are 
thus  attributed  to  the  presence  of  the  torsion  factor  T0=T/k.  Proceeding 
next  to  rewrite  Equation  (7)  completely  in  terms  of  the  surface  ray  coordi- 
nates, one  obtains 


*(«>♦!?  + 


• 'n  Tq  jL-  [uU)-v(!-)])d  G(ks) 


As  in  Equation  (27),  it  is  noted  that  D in  Equation  (29)  is  unity  for  the 
cylinder.  When  the  torsion  T Is  allowed  to  vanish,  the  above  result  in 
Equation  (29)  for  the  cylinder  also  reduces  exactly  to  the  corresponding 
result  in  Equation  (15)  for  the  sphere.  One  may  next  rewrite  Equation  (18) 

In  the  same  way  that  Equation  (7)  was  rewritten  to  obtain  (29).  Thus 

+ n'tT0  jL  [u(5)-vU)])d  G(ks).  (30) 

Finally,  Equation  (20a)  may  be  rewritten  as 

Ej(pNMYor2Pe-n'n({vU)-  ^ vU)+(jL-)2uU)J +Tq  ^ [u(c )-v(c)]|  D G(ks), 
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after  noting  that  f-l=T  for  the  circular  cylinder.  Again  the  factor  D In 

Equations  (30)  and  (31)  Is  unity  for  the  circular  cylinder  and  It  Is  Introduced 

symbolically  in  the  above  equations  to  facilitate  a direct  comparison  between 

the  cylinder  and  the  sphere  solutions.  Clearly,  Equation  (30)  reduces  directly 

to  the  sphere  result  of  Equation  (22)  when  the  torsion  is  allowed  to  vanish;  on 

the  other  hand  Equation  (31)  reduces  exactly  to  the  sphere  result  of  Equation 

(24)  if  T (and  hence  TQ)  is  set  equal  to  zero  except  for  the  term  (^-)^uU)G(ks) 

which  is  different  from  the  term  (j^r)2v(?)G(ks)  present  in  the  sphere  result. 

However,  the  dominant  differences  between  the  sphere  and  the  cylinder  solutions 

are  clearly  attributed  to  the  presence  of  torsion  associated  with  the  surface 

rays.  The  only  other  differences  that  exist  between  the  cylinder  and  sphere 

solutions  involve  terms  with  different  Fock  functions;  i.e.,  they  involve  terms 

of  the  type  (j^r)2u(s)G(ks)  and  (^-)2v(?)G(ks) . One  of  these  terms  occurs  in 

the  sphere  solution;  whereas,  only  the  other  occurs  in  the  cylinder  solution, 

and  vice  versa.  Notice  that  both,  (^-)2u(5)G(ks)  and  ( jg-)2v(?)G(ks)  behave 

essentially  as  (|j)2G(ks)  when  ks  is  small;  furthermore,  these  terms  of  order 

(^-)3  are  important  only  for  very  small  ks  values.  Consequently,  the  failure 

of  the  cylinder  solution  to  exactly  reduce  to  the  sphere  solution  to  terms  in 
1 3 

%>  • when  T is  set  to  zero,  is  not  a serious  problem  from  the  point  of  view 
of  being  able  to  generalize  these  canonical  solutions  to  the  arbitrary  convex 
surface.  The  following  heuristic  procedure  is  proposed  for  generalizing  the 
cylinder  results  in  Equations  (27),  (29),  (30),  and  (31),  and  the  corresponding 
sphere  results  in  Equations  (11),  (15),  (22),  (24)  to  treat  the  arbitrary 
convex  surface.  The  results  in  Equations  (27),  (29),  (30),  and  (31)  for  the 
cylinder  may  be  assumed  to  be  applicable  for  treating  an  arbitrary  convex 
surface  as  well  via  the  local  properties  of  propagation  and  diffraction  at 
high  frequencies, provided  the  geometrical  properties  of  the  arbitrary  surface 
are  slowly  varying  with  respect  to  the  wavelength.  The  geodesic  ray  paths 
on  an  arbitrary  convex  surface  in  general  possess  non-zero  torsion.  The 
effects  of  torsion  on  the  surface  ray  field  are  then  accounted  for  by  the 
factor  Tq  present  in  Equations  (27),  (29),  (30),  and  (31).  The  torsion  T 
and  curvature  < appearing  in  TQ  of  Equation  (28a)  are  given  by  Equations 
(28b)  and  (28c),  respectively  for  a circular  cylinder.  For  an  arbitrary 
convex  surface  one  may  write 


27 


(32a;32b;32c) 


sin2«5/ 


2 2 

k = cos  6<^+s1n  Sicg 


with 


£-t.|  = cos6  ; fc-Tg  = sin<s 


(32d;32e) 


In  which  and  t2  are  now  the  principal  directions  at  any  point  on  an 
arbitrary  convex  surface,  and  k ^ and  are  the  Principal  surface  curvatures. 
It  will  be  assumed  for  the  sake  of  definiteness  that 


K1  - k2 


<* 

Ki  >0 
1 

K«  > 0 

. C 


(32f) 


The  result  for  torsion  (or  T)  in  Equation  (32b)  which  applies  to  a geodesic 
surface  ray  path  on  an  arbitrary  convex  surface  has  been  obtained  in  a simple 
manner  through  the  use  of  differential  geometry;  whereas,  the  result  in 
Equation  (32c)  for  the  curvature  (or  K=l/pg)  is  simply  a statement  of  Euler's 
theorem.  For  surfaces  with  slowly  varying  geometrical  properties,  one  may 
symmetrically  split  the  factor  TQ  between  the  launch  (or  source)  point  and 
the  receiving  (or  field)  point  on  the  surface  to  represent  an  averaging 
type  effect.  Such  a symmetrical  split  is  also  essential  for  preserving 
reciprocity.  Thus* 
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tTpUT 


w 


OTTO 


(33) 


Next,  one  introduces  the  dimensionless  shape  factors  Ag  and  ac  to  appropriately 
weight  the  (^)2u(c)  6(ks),  and  the  -f^-)2v(5)S(ks)  type  terms,  such  that  only 
one  or  the  other  type  term  is  correctly  present  when  the  results  for  the 
arbitrary  convex  surface  are  specialized  to  either  the  circular  cylinder  or 
the  sphere.  In  particular,  this  would  require  that  A$  and  Ac  have  the  follow- 
ing properties: 


* The  - sign  is  chosen  in  (33)  if  - is  negative  at  P^  and/or  PN;  otherwise 
the  + sign  is  chosen 


i 


28 


(34a) 


A.+A  = 1 
s c 


a$M 


1 , sphere 
0,  cylinder 


(34b) 


fl,  cylinder 


|0,  sphere 


(34c) 


As  will  be  seen  later,  the  property  in  Equation  (34a)  is  also  essential  for 
recovering  the  planar  surface  result  as  k-*0.  A simple  choice  of  A$  would 
be  to  define  it  as  the  ratio  of  the  principal  curvatures;  thus. 


K1 

A = — = <> 
s <2 


1,  for  a sphere  (since  <^=<2) 
^0,  for  a cylinder  (since  k-j=0) 

From  Equation  (34a),  ac  must  then  be  defined  as 


(34d) 


K0-K 


2 1 _ 


1,  for  a cylinder  (since  iq=0) 
0,  for  a sphere  (since  k|=k2) 


(34e) 


Furthermore,  the  above  A and  A will  also  be  symmetrically  split  between 


the  launching  and  receiving  points  to  preserve  reciprocity  and  to  represent 
an  averaging  type  effect  for  surfaces  with  slowly  varying  geometrical  pro- 
perties. It  follows  that. 
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KlfW 

k^TpJ|T 


(PJ 


1 ' N 


A =1 -A  (with  A as  in  Equation  (35a) 

v 5 b 


(35a;35b) 


The  shape  factors  A$  and  ^c  essentially  serve  to  smoothly  interpolate  between 
the  cylinder  and  sphere  solutions;  this  role  of  the  shape  factors  together 
with  the  relation  in  Equation  (35a)  guarantees  that  the  conditions  set  forth 
in  part  (a)  are  satisfied  by  the  solution  for  the  arbitrary  convex  surface. 
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As  is  evident  from  the  cylinder  and  the  sphere  solutions,  the  variation 
of  the  surface  field  between  the  launching  and  receiving  points  is  primarily 
governed  by  the  Fock  type  functions  u{$)  and  v(£).  In  the  case  of  the  cylinder, 

the  Fock  parameter  s = — = ms<;  whereas,  s = ~ for  the  sphere.  This  Fock 

pg 

parameter  £ may  be  readily  generalized  in  the  usual  manner  to  treat  the 
arbitrary  convex  surface  case  by  employing  the  local  properties  of  wave  pro- 
pagation at  high  frequencies  as  follows. 


in  which  the  integral  is  evaluated  over  the  geodesic  surface  ray  path  from 
PN  PN*  **  1S  assume<*  development  of  the  result  in  Equation  (36) 

that  m and  hence  vary  slowly  with  a change  in  the  geodesic  arc  length; 
i.e.,  the  surface  properties  are  slowly  varying.  Besides  the  dependence  on 
pg  (or  1/ic),  the  variation  of  the  surface  field  between  and  P^  also  depends 

2 2-1 

upon  the  surface  radius  of  curvature,  p^^sin  6+k2cos  s)~  in  the  bi- 
normal direction  to  the  surface  ray  path;  this  dependence  occurs  through  the 
presence  of  the  spatial  factor,  D which  is  related  to  the  surface  ray  divergence 
factor.  For  an  arbitrary  convex  surface,  it  can  be  shown  that 


where  dij>0  is  the  angle  extended  by  the  surface  ray  tube  at  the  launching 
point  P^,  and  likewise,  dip  is  the  angle  extended  by  the  same  ray  tube 
(or  pencil)  at  the  receiving  point  P^  as  shown  in  Figure  4.  Also 
illustrated  in  Figure  4 is  the  caustic  distance,  pc  of  the  wavefront 
associated  with  the  surface  ray  tube.  The  caustic  distance  pc  is  the 
geodesic  (or  tangent)  radius  of  curvature  of  the  geodesic  circle  at  P^. 
It  can  be  shown  through  the  use  of  differential  geometry  that. 


where  E,  F and  G are  the  usual  "first  fundamental  coefficients"  which 
arise  in  the  development  of  the  first  fundamental  form  associated  with 
the  differential  arc  length  along  a curve  on  a surface.  It  is  easily 


verified  that  p =s  on  a cylinder  and  also  on  a plane;  whereas  p=atane-atan( 

C 

for  a sphere  (see  Figure  3).  This  concludes  the  discussion  corresponding 
to  the  assumptions  and  postulates  of  part  (b). 


p * CAUSTIC  DISTANCE 
C FROM  Pc  TO  PN 


Figure  4.  Caustic  distance  associated  wit  the  spreading 
of  the  surface  ray  field. 
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A.  Surface  Fields  of  an  Infinitesimal  Magnetic 
Current  Source  on  a Convex  Surface 

In  this  section,  an  asymptotic  solution  is  presented  for  the  field  on 
an  arbitrary  convex  surface  excited  by  an  infinitesimal  magnetic  current  source. 
The  construction  of  this  solution  is  heuristic  in  nature,  and  it  follows 
directly  via  the  postulates  and  assumptions  outlined  in  parts  (a)  and  (b)  above. 
An  analogous  solution  for  the  field  on  an  arbitrary  convex  surface  excited  by 
an  infinitesimal  electric  current  source  will  be  presented  later  in  Section  B. 
Thus,  the  field,  for  an  arbitrary  convex  surface  becomes  via  Equations  (11), 
(27),  and  the  postualtes  in  parts  (a)  and  (b),  the  following. 


where  the  dyadic  ~y  is  given  by 


(39a) 


% • (^jp-  feMO+D2(k)2[AsuU)+Acv({)K  |jCu(5)-v(s)]] 

+ Vt  (°2  is  v(?)+  is  u(c)"2(ii')  l\u(5)+Acv(0]}  + 

+ (K+b't)  jro  *L  [uU)-vU)]j)  D 6(ks)  . 


(39b) 


The  quantities  £,  D,  A , a * T and  G(ks)  are  as  defined  in  Equations  (36), 
(37),  (35a),  (35b),  (33),  and  (4),  respectively.  One  similarly  obtains  E^J 
for  an  arbitrary  convex  surface  as 


* Pm‘^i 


(40a) 


where  the  dyadic  Lm  is  given  by 


* (y0)"'(b'n  f[l-  jL.]v(s)+T2  jL  [u(0-*(5)]}  * 


+ {To  fe  Cuu>-vun})D  s(ks)  . (w) 

It  is  easily  verified  that  Equation  (39)  reduces  to  Equation  (27)  for  a 
circular  cylinder  when  and  Kg-l/a  (a-radlus),  and  It  reduces  to  Equation 
(11)  for  a sphere  when  Kj-tcg^l/a  To*0.  Likewise  Equation  (40)  reduces  to 
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Equation  (29)  when  it  is  specialized  to  the  circular  cylinder  case,  and  it 
reduces  to  Equation  (15)  when  it  is  specialized  to  the  sphere  case.  Further- 
more, both  Equations  (39)  and  (40)  properly  reduce  to  the  corresponding  results 
for  a planar  surface  when  <->0. 


Surface  Fields  of  an  Infinitesimal  Electric  Current 
Source  on  a Convex  Surface 


A 

The  field  Ht  for  an  arbitrary  convex  surface  becomes  via  Equations  (22) 
and  (30),  and  the  postulates  in  parts  (a)  and  (b),  the  following. 


"HV ' vre  • <4,a) 

where  the  dyadic  Te  is  given  by 

V - <V''("'6{n-  fcr  [»(«)-*(«)]}  ♦ 

S't  |t0  jL  [uUl-vU)]})  D G(ks)  . (41b) 

Likewise,  the  field  E^  for  an  arbitrary  convex  surface  is 

E^(Pn)  - pe-Ze  , (42a) 

where  the  dyadic  7g  is  given  by 

Ze  = - (Y0)“2  n'n^CvU)-  ^-v(5)]+(^-)  [asvU)+Acu(0]  + 

+ T2  jL-  [u(c)-v(0])  D G(ks)  . (42b) 

Again,  it  is  easily  verified  that  Equation  (41)  reduces  to  Equation  (30)  when 

it  is  specialized  to  the  circular  cylinder  case,  and  it  reduces  to  Equation 

(22)  when  it  is  specialized  to  the  sphere  case.  Also,  Equation  (42)  reduces 

to  Equation  (31)  when  specialized  to  the  circular  cylinder  case,  and  it  reduces 

to  Equation  (24)  when  specialized  to  the  sphere  case.  Finally,  both  Equations 

(41)  and  (42)  reduce  to  the  corresponding  results  for  a planar  surface  as  k-*0. 

It  Is  noted  that  for  a convex  surface  of  revolution  in  which  the  source 

p”  or  p„  on  the  surface  is  located  on  the  axis  of  revolution,  there  exists  a 
re  m 


i 
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caustic  of  the  surface  rays  at  a point  which  is  antipodal  to  the  source  point. 

In  this  case,  the  caustic  matching  functions  and  must  be  intro- 

duced in  the  results  of  Equations  (39),  (40),  (41),  and  (42),  in  exactly  the 
same  fashion  as  done  in  Sections  II-B  and  III-B  for  the  canonical  sphere  geo- 
metry. Also,  in  the  case  that  the  field  point  is  raised  from  on  the  surface 
to  the  point  P which  is  a small  distance  above  the  surface,  the  functions  u(?) 
and  vU)  in  Equations  (39),  (40),  (41),  and  (42)  must  be  replaced  by  FsU,0,y2) 
and  Fh(5,0,y2),  respectively.  In  addition,  if  the  source  Pg  at  P^  on  the  surface 
is  also  raised  slightly  above  the  surface  to  the  point  P',  then  Fs(c,y1,y2)  and 
FhU,y, ,y2)  must  be  used  instead,  as  indicated  previously. 

Finally,  if  the  surface  changes  continuously  but  not  sufficiently  slowly 
from  an  almost  spherical  surface  at  one  point  on  a geodesic  surface  ray  path  to 
an  almost  cylindrical  surface  at  another  point  on  that  geodesic  path,  then  the 
results  given  in  Equations  (39),  (40),  (41),  and  (42)  may  not  be  valid  because 
the  GTD  procedure  upon  which  they  are  based  is  not  strictly  valid  for  rapidly 
varying  surface  properties.  However,  if  the  change  in  surface  properties  is 
not  too  rapid,  then  it  would  be  interesting  to  study  if  one  could  heuristically 
take  such  moderately  rapid  changes  into  account  by  allowing  the  factors  Tq,  a$. 


and  ac  to  also  change  continuously  along  the  geodesic  path.  Thus,  one  may 


conjecture  that 


ro  s i Ip"  (™dS  ; As  = s L" 


(43a;43b) 


and 


Ac  = s 


) 

(1-(ki/k2^<J's 


(43c) 


It  is  noted  that  the  results  in  Equations  (39),  (40),  (41),  and  (42)  appear 
to  be  accurate,  but  they  still  are  being  tested  carefully;  as  such,  they  may 
be  subject  to  further  improvements  or  modifications  in  the  future.  Equations 
(43a;b;c)  are  added  simply  as  something  of  interest  to  study  in  the  future, 
and  it  is  not  recommended  that  they  be  used  until  sufficient  tests  on  their 
validity  have  been  completed. 


V.  NUMERICAL  RESULTS  AND  DISCUSSION 


In  this  section,  general  expressions  will  be  indicated  for  the  mutual 
coupling  between  a pair  of  antennas  located  on  a smooth,  perfectly-conducting 
convex  surface  of  any  shape.  These  expressions  will  of  course  make  use 
of  the  results  developed  in  Section  IV.  The  type  of  antennas  considered  in 
this  paper  are  thin,  rectangular  slots  which  are  less  than  half  a wavelength 
long,  and  thin  stubs  which  are  less  than  a quarter  wavelength  long.  The 
reason  for  selecting  such  short  and  thin  slot,  or  stub  antennas  is  that 
the  shape  of  the  field  distribution  in  the  slot  or  the  current  distribution 
along  the  stub  may  then  be  assumed  to  be  that  of  the  dominant  mode. 
Furthermore,  a short  circuit  placed  in  a waveguide  feeding  the  slot  such 
that  it  is  an  integral  number  of  "dominant  mode"  half  wavelengths  from  the 
slot  aperture  would  effectively  produce  a short  circuit  at  the  aperture. 

Thus,  a surface  field  which  propagates  along  a geodesic  surface  ray  path 
over  such  a short  circuited  aperture  would  for  all  practical  purposes  be 
unaltered  by  the  presence  of  the  slot  since  the  scattering  of  the  surface 
field  by  the  short  circuited  slot  would  then  be  vanishingly  small.  This 
fact  is  particularly  useful  if  one  is  interested  in  the  calculation  of  the 
mutual  admittance  between  a pair  of  slots  in  a slot  array  environment.  A 
similar  simplification  results  in  the  calculation  of  the  mutual  impedance 
| between  a pair  of  thin  and  short  stub  (or  monopole)  antennas  in  a stub 

| array  environment  since  such  an  open  circuited  stub,  placed  in  the  path  of 

tr 

the  surface  field  would  perturb  this  field  negligibly.  It  is  noted  that 
the  conditions  imposed  by  the  above  assumptions  on  the  size  of  the  radiating 
elements  are  not  stringent;  in  fact,  these  conditions  are  commonly  met  in 
practice.  The  calculation  of  the  mutual  coupling  between  a pair  of  slot  or 
wire  type  radiating  elements  on  a convex  surface  which  do  not  satisfy  the 
above  assumptions  and  conditions  can  still  be  performed  on  the  basis  of 
the  results  developed  in  Section  IV  provided  the  aperture  field  distribution 
(for  the  case  when  the  radiating  element  is  a slot  antenna  which  may  not  be 
sufficiently  small)  and  the  current  distribution  along  a wire  (for  the  case 
of  a stub  antenna  which  may  not  be  sufficiently  short  and  thin)  is  known  in 
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the  coupled  environment.  In  the  latter  cases,  one  could  employ  measurements, 
or  numerical  techniques  (such  as  the  method  of  moments)  for  solving  the 
coupled  integral  equations  for  these  problems  in  order  to  obtain  the  field 
or  current  distributions  on  these  radiating  elements.  These  integral  equations 
could  be  simplified  so  that  the  only  unknowns  are  the  fields  on  the  antennas 
and  not  on  the  electrically  large  structure  upon  which  the  antennas  are 
located,  by  appropriately  employing  the  surface  field  dyadics  introduced  in 
Section  IV,  in  the  construction  of  these  integral  equations.  These  dyadics 
would  actually  constitute  asymptotic  approximations  to  the  appropriate 
surface  dyadic  Green's  functions  or  their  spatial  derivatives  which  con- 
stitute the  kernels  of  these  integral  equations. 


Some  numerical  results  will  be  presented  in  this  section  for  the 
calculation  of  the  mutual  admittance  between  a pair  of  rectangular  slots  in 
a perfectly-conducting  circular  cylinder;  both,  axial  and  circumferential 
slots  will  be  considered.  These  numerical  results  which  are  based  on  the 
development  in  Section  IV,  will  be  compared  with  the  corresponding  numerical 
results  obtained  by  others  [11,12]. 


The  general  expressions  for  the  mutual  coupling  between  a pair  of 
radiating  elements  on  a smooth,  perfectly-conducting  convex  surface  are 
presented  in  the  following  discussion.  In  particular,  let  antennas 
designated  (l)  and  (2)  constitute  a pair  of  radiating  elements  on  a convex 
surface.  The  mutual  coupling  between  these  antennas  may  be  calculated  via 
an  application  of  the  reaction  theorem  [2].  Thus,  if  antenna  (2)  is  an 
aperture  (or  slot),  the  dominant  mode  current  ^ induced  in  this  aperture 
by  antenna  (l)  when  antenna  (D  is  short  circuited  is  given  by  [2] 


Ifds2  H21  *^2 

s2 


v22 


(44) 


,! 
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where  ^2  = ^2  x n2’  and  ^2  1S  dom^nant  mode  electric  field  distribution 
in  the  aperture  of  the  slot  antenna  (?)  with  ng  being  the  unit  outward 
normal  vector  to  the  surface  of  the  aperture  of  antenna  (§).  The  inte- 
gration in  Equation  (44)  is  over  the  surface  area  S£  of  the  aperture 
corresponding  to  slot  antenna  d).  Furthermore, 


?e  = v22uy  Jar  cos<r  ux>  • 


In  Equation  (45)  above,  is  the  dominant  mode  voltage  associated  with 

the  dominant  vector  mode  function  e”  for  a rectangular  slot.  The  local, 

orthogonal  coordinates  u and  u whose  origin  is  at  the  center  of  the  slot 

* y 

are  indicated  in  Figure  5 together  with  the  slot  dimensions  d and  *, 


SLOT  ANTENNA 


luxl<  V2 

SLOT  DIMENSIONS:  <luyl  < d/2 

u.  *0 


* _ A w A 

uz  ■ uK*uy 


Figure  5.  Local  coordinates  associated  with 
the  slot  antenna. 

respectively.  The  field  is  the  magnetic  field  at  antenna  (f),  which 
is  produced  by  antenna  (T)  with  antenna  (2)  short  circuited.  From  the 
development  in  Section  IV,  it  is  easily  seen  that  when  antenna  (1)  is 


also  a slot,  then  *21  = II  dSi  where  the  integration  is  over  the 

S1 

surface  s-j  of  the  aperture  corresponding  to  slot  antenna  (l) , with  t/m  as 
given  in  Equation  (39b),  and  Fj  = ^ x n-j.  Of  course,  n-j  is  the  outward 
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unit  normal  vector  to  the  surface  area  of  the 
to  slot  antenna  0,  and 

ri = vn  • • 

v-j -j  is  the  dominant  mode  voltage  associated  with 
other  hand,  if  antenna  0 is  a stub,  then  F2-|  = 


aperture  corresponding 


(46) 


slot  antenna 

JJjTdv,  j,-t 


0 . On  the 
e,  where  the 


integration  is  on  volume  v1  just  encapsulating  stub  0,  with  as  given  in 
Equation  (41b),  and 


J1  = !11^=  Iiiuz6(ux)6(uy)1(u2^ 


(47a) 


In  Equation  (47a)  above,  1^  is  the  mode  current  on  the  stub,  which  is 

associated  with  the  vector  mode  function  F for  the  stub.  The  vector  F 

involves  the  2-D  Dirac  delta  function  S (u  )S(u  ) and  a distribution  I (u  ) 

x y z 

which  is  generally  assumed  to  be 
sin[k(*-|u  I)] 

I(uz>  = sin(kt)  ! 0 < u2  < t . (47b) 


The  local  orthogonal  coordinates  (u  ,u  ,u  ) whose  origin  is  at  the  base 

a y 4 

(or  feed  point)  of  the  stub  are  illustrated  in  Figure  6.  From  the  above 
discussion,  it  follows  that 
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, if  antenna 
0 is  also 
a slot, 


and 


(48) 
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STUB  ANTENNA 
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J JJ  dv2  *21  ‘7Z 


where  v2  is  the  volume  just  encapsulating  stub  antenna  0,  and  J2  * I22  F 
on  this  stub.  is  the  electric  field  at  stub  antenna  0,  which  is 
produced  by  antenna  0 when  antenna  0 is  open  circuited.  It  follows 

from  the  development  in  Section  IV  that  F21  = JJ  ds-j  if  antenna  0 

is  a slot,  and  r21  = JJT  dv-|  if  antenna  0 is  also  a stub.  The 
= = V1 

dyadics  Lm  and  Zg  are  given  in  Equations  (40b)  and  (42b).  Summarizing  the 
results  for  the  open  circuit  voltage  v^  at  stub  antenna  0,  one  has 


Ilf  dv2  JJ  dsi 


Stub  0 
Open  Circuited, 

O'"  l2i=0 


, if  antenna  0 
is  a slot. 


J/Jdv2  JJf  dvl  J1 


VT2 


21  Stub  0 

Open  Circuited, 


, if  antenna  0 
is  also  a stub. 
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It  is  evident  from  Equation  (53),  that  the  mutual  impedance  Z2-|  between 
a pair  of  stubs  (monopoles)  is  given  by  V^/I^  with  I2-,=0,  or 


» 


Z21  = Z12* 


(54) 
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Some  numerical  results  are  indicated  below  for  between  a pair 
of  identical  rectangular  slots  in  a perfectly-conducting  circular  cylinder. 
In  particular.  Figures  7 and  8 indicate  the  surface  field  on  a circular 
cylinder  due  to  an  infinitesimal  circumferential  slot;  whereas  Figures  9 
and  10  indicate  the  isolation,  s12  (related  to  Y12)  between  a pair  of 
axial  slots.  These  numerical  results  are  marked  as  OSU  on  the  plots. 
Furthermore,  these  results  are  compared  with  the  exact  (modal)  results 
as  well  as  those  in  [11,12];  in  particular,  the  results  of  [11]  are  marked 
as  PINY  in  these  figures;  whereas,  those  of  [12]  are  marked  as  UI  results. 
Additional  comparison  between  the  OSU,  PINY,  and  UI  results  for  Y-|2  are 
indicated  in  Tables  I and  II  for  the  circumferential  slots.  It  is  seen 
from  these  comparisons  that  firstly  the  OSU  results  agree  very  well  with 
the  results  based  on  the  exact  (modal)  solution  given  in  [12],  which  also 
contains  the  Hughes  modal  solution.  Secondly,  it  appears  that  the  results 
based  on  all  three  asymptotic  solutions;  namely  those  based  on  the  OSU, 

PINY  and  UI  solutions  are  for  all  practical  purposes  quite  accurate  and  not 
noticeably  different  except  when  ks  becomes  small,  and/or  the  pair  of  slots 
are  in  each  other's  paraxial  regions  for  the  case  of  circumferential  type 
slots.  Thus,  when  ks  becomes  very  small,  the  PINY  solution  becomes 

inaccurate  because  only  terms  up  to  and  including  — — appear  in  that 

(ksr 

solution;  whereas,  the  OSU  and  UI  results  which  include  the  effects  of 
1 3 

(j~)  terms  remain  accurate  even  for  ks  very  small.  The  OSU  and  PINY 
asymptotic  solutions  for  the  circular  cylinder  are  not  strictly  valid  within 
paraxial  regions  since  the  approximation  for  the  field  in  terms  of  Fock 
type  Airy  functions  cannot  be  completely  justified  in  these  regions.  On  the 
other  hand,  the  UI  solution  bypasses  this  difficulty  by  conjecturing  an 
asymptotic  solution  for  the  circular  cylinder  via  a heuristic  modification 
of  an  asymptotic  solution  for  the  sphere  geometry.  As  mentioned  previously 
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TABLE  I 


Y]2  Between  a Pair  of  Circumferential  Slots  in  a Circular 
Cylinder  of  Radius,  a=1.991";  frequency=9  GHz,  and  6=0  (E-plane) 


1 

Modal  Sol 

utions 

Asymptotic  Solutions 

1 

(Inches) 

UI 

UI 

OSU 

PINY 

Planar 

a*«* 

1 

.5 

- 7.27  dB 

1 9 ' 

- 7.31 

- 7.86 

- 6.46 

- 8.16 

r 

-72° 

-77.77° 

-69.71° 

-68.17° 

-66.85° 

t 

2 

-16.52 

-16.36 

-17.56 

-15.66 

mmm 

-117° 

-117° 

-115.67° 

-110.08° 

-117.88° 

-105.84° 

5 

8 

-26.95 

-26.49 

-26.54 

-28.78 

-25.51 

-28.97 

33° 

34° 

36.77° 

46.20° 

33.60° 

53.6° 

| 

i 

16 

-31.13 

-4° 

-31.31 

-.9° 

-34.28 

10.19° 

-30.04 

-4.2° 

-35.98 

19.96° 

40 

-36.60 

-115° 

-37.17 

-109.44° 

-41.29 

-96.97° 

-35.58 

-112.59° 

-43.93 

-83.17° 

i 

{ 

i 


d 


A *11 

= 0.9 

^ -H 

* 0.4 


I 

! 

3 

x 

i 


I 

i 
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TABLE  II 

Y-J2  Between  a Pair  of  Circumferential  Slots  in  a Circular 
Cylinder  of  Radius  a=1.991";  Frequency=9  GHz  and 


for  6=u/2;  zQ=0  (H-Plane) 


Modal 


-25.98  dB 
-77° 

-34.52 

108° 

-40.96 

58° 

-46.62 

-49° 


-34.63 

169.58° 

-41.32 

59.88° 

-47.08 

-47.85° 


totic  Solutions 

OSU 

PINY 

-26.07 

-75.73° 

-34.67 

170.07 

-41.76 

60.33° 

-46.92 

-47.55° 

-77.93 

-64.81° 

-35.72 

179.40 

-42.03 

68.31° 

-47.55 

-40.48° 

/ = 0.9" 
d = 0.4" 


PHASE  ( DEGREES  I 


20  L06l0  |SI2  I (dB) 


Figure  9.  Isolation  of  axial  slots  on  a conducting  cylinder 
a = 1.991";  Z0  = 1.50";  Frequency  = 9 GHz. 


20  LOG  iq  | S |2  ( (dB) 


4>  (AZIMUTHAL  SLOT  SEPARATION  IN  DEGREES) 


Figure  10.  Isolation  of  axial  slots  on  a conducting  cylinder, 
a = 1.991";  Z0  = 0;  Frequency  = 9 GHz. 
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(Section  I),  this  modification  of  the  sphere  solution  in  [12]  results  from  the 
introduction  of  an  extra  term  via  an  ad  hoc  procedure.  This  extra  term  improves 
the  accuracy  of  the  modified  sphere  result  when  it  is  specialized  to  the  circular 
cylinder  case;  and  the  improvement  is  especially  noticeable  in  the  paraxial 
regions  of  the  cylinder.  As  a result,  the  UI  solution  appears  to  be  quite  ac- 
curate even  in  the  paraxial  regions,  i.e. , in  the  immediate  vicinity  of  the 
axial  direction  on  a circular  cylinder  for  the  circumferential  slot  case. 

Although  the  OSU  and  PINY  solutions  for  the  circular  cylinder  are  not  strictly 
valid  within  te  paraxial  regions,  they  still  provide  accurate  results  for  the 
axial  slot  case,  and  they  are  reasonably  accurate  even  for  the  circumferential 
slot  case;  thus,  these  solutions  may  be  employed  even  within  the  paraxial  regions. 
Actually,  the  OSU  results  for  the  paraxial  regions  corresponding  to  the  circum- 
ferential slot  case  differs  from  the  exact  (modal)  results  by  less  than  10%  in 
amplitude,  and  by  less  than  5%  in  phase.* 


It  is  evident  from  the  results  in  Sections  II  and  III,  that  the  asymp- 
totic solutions  for  the  cylinder  are  somewhat  different  from  those  for  the 
sphere  as  one  might  expect,  since  the  two  geometries  are  different.  Also,  it 
is  noted  that  the  OSU  solution  for  the  circular  cylinder  case  (as  presented  in 
Section  II-A)  is  actually  identical  to  the  corresponding  PINY  solution  to  terms 

in  and  — — except  for  the  cross  term  (i.e.,  the  type  term)  whose 
ks  d d 

functional  form  is  different  in  the  OSU  and  PINY  solutions.  A cross  term  of 

this  type  is  absent  in  the  UI  solution.  Furthermore,  it  is  this  cross  term 

which  dominates  the  solution  in  the  paraxial  regions  for  the  circumferential 

slot  case.  Thus,  the  numerical  results  based  on  the  OSU  and  PINY  solutions 

are  slightly  different  in  the  paraxial  regions  for  the  circumferential  slot 

case.  While  the  UI  solution  surprisingly  yields  good  accuracy  for  the  circular 

cylinder  case,  it  raises  some  interesting  questions  which  were  also  indicated 

previously  in  Section  I.  In  particular,  the  UI  solution  when  specialized 

to  the  circular  cylinder  geometry  excited  by  an  axial  slot  will  contain 

both,  the  soft  and  hard  type  Fock  functions  u(?)  and  v(c),  respectively. 


♦It  is  expected  that  further  improvements  in  the  paraxial  region  field  calcu- 
lations are  possible  if  higher  order  terms  in  are  included  in  the  analysis. 

kpg  1 


However,  as  indicated  above,  the  leading  (or  lowest  order)  terms  in  ^ — which 


are  the  only  ones  that  have  been  retained  in  this  analysis  are  sufficiently 
accurate  even  for  paraxial  region  calculations. 
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However,  it  is  well  known,  that  only  the  potential  corresponding  to  the 
"hard"  or  v(c)  type  function  is  required  to  completely  describe  the  fields 
of  an  axial  slot  in  a circular  cylinder.  Secondly,  the  "ad  hoc"  procedure 
employed  to  modify  the  asymptotic  solution  for  the  sphere  in  order  to 
obtain  an  asymptotic  solution  for  the  circular  cylinder  geometry  may  need 
further  justification,  In  this  sense,  the  present  approximate  (OSU ) 
asymptotic  solution  for  the  circular  cylinder  is  not  based  on  a heuristic 
procedure  as  the  one  in  [12],*  and  it  may  be  viewed  as  an  extension  of  the 
approximate  PINY  asymptotic  solution  [11]  and  the  previous  OSU  solution  [9]. 

It  is  noted  that  the  OSU  asymptotic  solution  for  the  circular  cylinder 
could  be  improved  even  further  in  the  paraxial  regions  for  the  circumferential 
slot  case  by  employing  a different  asymptotic  evaluation  which  would  be 

more  accurate  in  these  regions;  however,  this  analysis  will  not  be  pursued 
at  the  present  time  since  the  solution  derived  in  this  paper  appears  to  be 
sufficiently  accurate  even  for  this  special  case. 


It  is  noted  that  In  addition  to  the  asymptotic  solutions  pertaining 
to  the  surface  fields  on  the  canonical  circular  cylinder  and  the  sphere 
geometries,  an  approximate  asymptotic  solution  is  also  heuristically 
constructed  in  this  paper  for  the  arbitrary  smooth  convex  surface.  Also, 
the  effect  of  torsion  associated  with  the  surface  rays  is  clearly 
identified  in  this  solution  through  the  presence  of  the  factor  T/k. 
Furthermore,  the  excitation  of  this  convex  surface  by  both  infinitesimal 
electric  and  magnetic  type  current  moments  are  considered  in  this  paper. 
Numerical  results  for  both  the  electric  and  magnetic  type  source 
excitation  of  spheroidal  geometries  will  hopefully  be  presented  in  the 
near  future  along  with  experimental  results  for  comparison  in  order 
to  test  the  validity  of  the  present  asymptotic  solution  for  the  arbitrary 
convex  surface. 


♦During  the  editing  of  the  final  version  of  this  paper,  we  have  received  from 
Professor  S.  W.  Lee  (at  University  of  Illinois)  a set  of  handwritten  notes  by 
Professor  J.  Boersma  which  deal  with  a new  and  more  rigorous  asymptotic 
solution  (that  is  valid  within  the  paraxial  regions)  for  the  surface  fields  of 
slots  on  cylinders.  These  notes  are  soon  expected  to  be  published  as  a UI 
report  [16].  It  is  interesting  to  note  that  the  solution  in  our  paper  agrees 
exactly  with  the  one  in  Professor  Boersma 's  notes  to  all  orders  in 

(rr)  and  to  leading  terms  in  (t~),  even  though  it  is  derived  via  a different 

Kp  Q 

asymptotic  procedure.  The  effect  of  torsion  associated  with  the  surface  rays 
has  not  been  specifically  identified  in  Prof.  Boersma' s notes. 
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APPENDIX  I 

ASYMPTOTIC  ANALYSIS  OF  THE  SURFACE  FIELDS  EXCITED  BY 
INFINITESIMAL  ELECTRIC  OR  MAGNETIC  CURRENT  SOURCES 
ON  A PERFECTLY-CONDUCTING  CIRCULAR  CYLINDER 


The'geometry  of  this  canonical  problem  is  illustrated  in  Figure  1. 

The  case  of  a tangential  magnetic  current  source  at  will  be  analyzed 
first;  the  analysis  for  a radial  electric  current  source  at  P^  will  be 
performed  subsequently. 

Let  a magnetic  current  source  R"  generate  the  surface  fields  Ejj  and 
^ on  the  surface,  where  is  the  radial  component  of  the  electric  field, 
and  is  the  tangential  component  of  the  magnetic  field.  The  source  R"  is 
defined  as 

■ Pm  «(  W - r'  |)  , (A— 1 ) 

where  F'  is  the  position  vector  of  and  F is  an  arbitrary  position 
vector.  The  quantity  pm  is  the  strength  of  the  infinitesimal  magnetic 
current  moment  which  is  oriented  tangential  to  the  surface. 

The  fields  E™  and  hT  can  be  constructed  from  a suitable  set  of  electric 
n t _ _ 

and  magnetic  vector  potentials  F and  A,  respectively  [14].  Following  the 
development  in  [14]  for  an  eJut  time  dependence  which  is  assumed  and  sup- 
pressed, one  obtains, 


where  (k^p)  is  the  Hankel  function  of  the  second  kind  and  of  order  n. 
The  kt  appearing  in  the  argument  of  this  Hankel  function  is  defined  by 

k2  = k2-h2  . ( 

The  unknown  coefficients  an(h)  and  bn(h)  in  the  above  expansions  may  be 
obtained  by  an  application  of  the  electromagnetic  boundary  conditions  on 
the  surface.  From  the  usual  relationships  between  the  vector  potentials 
and  the  electromagnetic  fields,  one  obtains  [14] 

H* = jk  ‘I? + k2)  F*  ( 


„ ..ii+  i 

<t>  ap  jkzQp  a<t>  3z 


E . 1 

p jkVQ  3p3z 


In  the  above  equations,  zQ  = free  space  impedance,  and  YQ  = l/zQ. 

Then,  + Hzz  , and  E^  = p Ep.  In  Equations  [A-5),  (A-6)  and 
(A-7),  the  fields  are  not  necessarily  on  the  surface.  For  convenience 
of  analysis,  p"  is  decomposed  as  follows. 


p = pa  t'  + pc  j' 

m Km  1 Hm  2 


T]  ! Z at 

*2  = ♦ at  P^ 


From  the  boundary  condition  Efnxn  = FTatp  = a,  where  E"1  is  the  total 
electric  field  and  is  as  in  Equation  (A-l),  it  can  be  shown  that 
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(A-10) 


In  order  to  evaluate  Equations  (A-5),  (A-6)  and  (A-7)  asymptotically  for 
large  kta,  one  begins  by  employing  the  usual  Watson  transformation  [5]  to 
Equations  (A-9)  and  (A-10).  Thus, 

«*.■  fc  f -*  T„ 

(A-ll) 

and 

F=zF  = L f dheJhz  f dvljf-y]  Hy.nktP>  H 

**  L Lj,  ir^]  ’ (kta)l*io^  > J 


- f dheihz  f ^ dv^j-Tl  H|;l!ktP>  [" | ^e-3'>*+e-j''(2’-♦>)e-J2’'l', 


(A-12) 


It  is  assumed  that  ka  is  sufficiently  large  so  that  one  may  retain  only 
the  t=o  term  in  Equations  (A-ll)  and  (A-12),  as  the  t>0  terms  correspond  to 
the  field  which  multiply  encircles  the  cylinder  and  contributes  insignifi- 
cantly for  large  ka.  The  contours  of  integration  in  the  complex  h and  v 
planes  are  shown  in  Figures  A-I,  and  A-II,  respectively.  The  dominant 
contribution  to  the  v integrals  in  Equations  (A-ll)  and  (A-12)  will  occur 
for  vM)(kta)  in  the  deep  shadow  region  and  also  in  the  shadow  boundary 
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transition  region.  Since  the  surface  fields  are  of  interest,  and  W lies 
on  the  cylinder,  the  field  point  is  always  in  the  shadow  region.  Thus, 
one  may  approximate  the  Hankel  functions  by  the  Fock  type  Airy  functions 
[6,15]  for  vM)(kta)  as  follows. 


Hv2><kta)  "2M 

Mr 

Hv2)(kta)  * ^ (mt2)  W2(t) 


(A-l 3) 


( A— 1 4) 


v = k^a  + m^r  ; 


k.a  1/3 

”t  * <-f  > 


(A-15a;A-15b) 


It  is  noted  that  W2(x)  is  the  derivative  of  W2(x)  with  respect  to  the 
argument  x.  The  functions  W^(t),  Wg(x)  and  their  derivatives  are  discussed 
in  Appendix  V.  Introducing  the  transformation  of  Equation  ( A-l 5a ) into 
Equations  (A-l 1 ) and  (A-l 2 ) together  with  the  approximations  in  Equations 
(A-l 3)  and  (A-l 4 ) , one  obtains. 
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(A-l 7) 


The  superscript  "+"  on  A*  and  F*  in  Equations  (A-l 6)  and  (A-l 7)  signifies 
that  only  the  term  corresponding  to  e"JV*  in  Equations  (A-l  1 ) and  (A-l 2 ) are 
considered  for  convenience;  the  expressions  for  the  e“jv(2w-<l>)  terms  are 
similar  and  may  be  designated  by  A^  and  F^  with  the  understanding  that 


Az  - Az  + Az  *. 


Fz  = Fz  + Fz  • 


(A-18a;A-18b) 


In  order  to  evaluate  the  integrals  with  respect  to  the  variable  h, 
it  is  convenient  to  introduce  the  usual  polar  transformations  given  by, 

h = ksina,  k*  = + kcosa;  d = a<f  = s cosi<»;  z = s sinip.  (A-19a;A-19b 
x A=19c;A-19d) 

Thus,  Equation  (A-16)  becomes  the  following  via  Equation  (A-19). 


jkY  D~ 
j^,  _ orm 


|p=a  4rr  a 


J da  e-jkscos(“-*!  / d,  -(A.20) 


The  contour  of  integration,  CQ  is  indicated  in  Figure  A- III.  The  dominant 

lm  a ! \ 


Re  a 


as  IS  THE  SADDLE  POINT 
LOCATION 


i • I 

i i 

I I 

Figure  A- III.  Contours  Ca  and  CSDp  in  the  complex  a plane. 

contribution  to  the  integral  over  a occurs  for  a * a$  = ip  for  large  ks  where 
a$  is  the  saddle  point  in  the  complex  a plane.  The  integral  over  t yields 
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a slowly  varying  function  of  a in  the  neighborhood  of  a$.  The  contour 
Ca  may  be  continuously  deformed  into  the  steepest  descent  path  C^p  since 
no  singularities  are  enclosed  (for  i/>  < tt/2) , and  the  contributions  from  the 
paths  connecting  Cq  with  Cs[Jp  are  vanishingly  small.  Therefore,  Equation 
(A-20)  becomes 


3 Al  jkY 

_JL  %. 

p=a  4it  a 


. OTT 

rr3/2^  J dcte-J'kscos^-^i 


(A-21 ) 


where  the  Fock  type  function  u(0  is  discussed  in  Appendix  V;  it  is  defined 


3/2  1 f 

■ * 5 ' k j. dT  wr 


. 


(A-22a) 


kp  1/3 


5 = 5 m = ( Z1 1 ; Pq  = 


(A-22b;A-22c;A-22d) 


cos  > 


The  geometrical  interpretation  of  the  angle  ij>  is  illustrated  in  Figure  1 
in  terms  of  the  angle  6 where 

6 = * * = " cos'1  (t'. . (A-23) 


It  can  be  shown  that  the  leading  term  in  the  steepest  descent  approximation 
for  large  ks  yields  the  following  relationship  for  the  CSDp  integral  in 
Equation  (A-21). 


J do  e"J"kscos(o",,')  <v  ^ 


(A-24) 


It  is  interesting  to  note  that  the  left  hand  side  of  Equation  (A-24)  is 
directly  related  to  another  integral  as  follows 
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1 

2J 


2jkcos^ 
\ ird 


j da  e-jkscos(“-*).  (A-25a) 
SDP 

for  |ktd|»T 

with  kt  = >k2-h2. 


Furthermore,  it  is  well  known  that  the  left  hand  side  of  Equation  (A-25a) 
is  exactly  given  by 
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2T 


oo 

I 


dh 


H‘Z)(ktd)eJhz  = 


,-jks 


■4 


? — 7 
<r+z* 


(A-25b) 


Now  the  approximation  arrived  at  on  the  right  hand  side  of  Equation  (A-24) 
is  identical  to  the  result  on  the  right  hand  side  of  Equation  (A-25b).  Thus, 
the  leading  term  of  the  asymptotic  approximation  to  the  CSDp  integral  in 
Equation  (A-25a)  is  identical  to  the  exact  result  in  Equation  (A-25b). 

Based  on  these  observations,  it  would  appear  that  the  final  result  in 
Equation  (A-24)  is  therefore  valid  even  for  ks  small.  Equation  (A-21)  may 
now  be  written  as, 

I 

* - pj  f ir  u(C)  G(ks)  (A-26) 

p=*a 

wlt’i 


G ( ks ) = 


e-jks 
~~W  ’ 


1/cos2^  * l/sin^6  . 


(A-27  a) 
(A-27b) 


Similarly,  one  can  show  that  Equation  (A-17)  becomes 
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• J da  e -jkscos  (a-ij>) 
CSDP  _ 


(A-28) 


where  the  Fock  type  function  v(0  is  discussed  in  Appendix  V;  it  is  also 
defined  below. 

HT  f W0(t) 


(77  r tMrj 

,u>  = Jte  J dT  snry 


(A-29) 


The  third  term  on  the  right  hand  side  of  Equation  (A-28)  is  of  higher  order 
in  {7  as  compared  to  the  other  two  terms;  hence,  this  term  will  be  discarded 
as  only  the  lowest  order  terms  in  are  retained  in  this  analysis.  Con- 
sequently, Equation  (A-28)  may  be  written  as 


F+  *!jn  ml  Ev(c)  f dae-jkscos(a-if.) 
Z p=a  4 ir  pg^c  / 


m kcots 


Jfv(5)  / dae-Jkscos‘“-ri  . CA— 30) 

/■ * 


From  Equation  (A-24),  and  Equation  (A-30),  one  obtains 
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|p*a 


-iks  Pc 

V--^C0t6  v(0^- 


(A-31) 


From  Equations  (A-5),  (A-6),  (A-26)  and  (A-31), 


. , /V  2\  * - 8ftj  i *'  K 

— twp=8 


(A-32) 


with  the  understanding  that  FI^(P^)  in  Equation  (A-32)  is  associated  only 
with  the  A*  and  F+  terms;  the  total  mus^  ^ obtained  from  (A*+A~) 

and  (F*+F“).  However,  the  form  of  Ti!£(PN)  associated  with  the  (A”,F~) 
terms  will  be  identical  to  that  associated  with  (A*,F*).  In  order  to 
simplify  the  derivatives  in  Equation  (A-32)  one  employs  the  following 
relationships. 


||-  = cosfi  = sin^;  14  = r-  cos6;  = sin6;  = ~ sinfi. 


32  p. 


(A-33a) ;(A-33b) 
(A-33c); (A-33d) 


Employing  Equation  (A-33)  into  (A-32),  and  retaining  only  terms  to  lowest 
order  in  one  obtains 

HJ(PN)  * Pm^^1[sin26+  j^l-  ^-)(2-3sin267j  v(?)G(ks) 

+ Pm*;2^l[2cos6sin6Y"  Is0"  fej]vU)G(ks)  + 

+ Pm-^t20in6cos6^-  §(1-  ^|jv(C)G(ks) 

+ Pm^2^|^s26+  fe*1"  j|-)(2-3cos26j]  v(0 

+ [u(c)-v(5)^G(ks)  . (A- 

The  above  result  has  the  interpretation  that  the  field  propagates  along 
the  helical  geodesic  surface  ray  path  from  to  PN  as  shown  in  Figure  1. 

A second  helical  geodesic  ray  path  (not  shown  in  Figure  1)  which  traverses 


(A-34) 
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the  cylinder  surface  in  a sense  opposite  to  that  shown  in  Figure  1 also 

contributes  to  the  total  field  at  P^;  the  latter  contribution  results  from 

the  (A",F")  terms»  and  it  is  identical  in  form  to  that  in  Equation  (A-34), 

except  that  s,  a and  e are  now  associated  with  this  second  path.  It  is 

noted  that  t,  and  t,  correspond  to  zL  and  z|D  , respectively;  whereas, 

1 N ,PN 

r'2  and  t2  correspond  to  $|p,  and  5>|p  , respectively. 

Next,  one  may  evaluate  Ep  due  to  FT  via  Equations  (A-7),  (A-26) 
and  (A-31 ) as  follows. 

ep*-  7^7) + ^ 


< 3 [(£>e'Jkrl 

■ 7T  vU)-rr 


(A— 35) 


where  d=a<j>.  One  may  simplify  Equation  (A-35)  by  employing  the  relations 
in  Equation  (A-33);  thus,  retaining  terms  to  lowest  order  in  1/ka, 


h zl-  ^ ;iV 

‘ (ksslnsrj  (ks)k  $in‘« 


■ e-jks  cos « 
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cos6  e-jks 


sin  6 


k2p®  P . -jks 
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in6  , 


(A-36a) 


or  in  the  vector  notation  of  Section  II, 
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r 


^n(PN)^(Yo)"lpm‘^isin6+^2cos6]"(1" 

+ Tgcosfin  jL{u(c)-vU)]  f la(ks). 


(A-36b) 


One  may  now  develop  similar  expressions  for  the  electromagnetic 
surface  fields  of  a radial  electric  current  I on  a perfectly-conducting 
circular  cylinder.  One  defines  T in  a manner  anologous  to  that  in  Equation 
(A— 1 ), 


Pe5(|r-r'|). 


(A-37) 


It  is  noted  that  pg  is  the  strength  of  the  infinitesimal  electric  current 
moment  which  is  oriented  normal  to  the  surface  at  F' . The  tangential  magnetic 
field  ^(PN)  due  to  T may  be  simply  obtained  by  employing  the  reciprocity 
theorem  to  the  fields  of  M and  T;  namely. 


nj-pJ 

t m 


at  p 


- E^p 
n He 


m 


(A-38) 


at  p. 


Consequently,  from  Equations  (A-36b)  and  (A-38),  one  obtains 
^(PN)'v-(Y0)"1pe‘^'[-T1sin5+;?2cos6](l-  jL-)vU) 

+n'cos6  t2  jL{uU)-v(C)]  f^G(ks).  (A-39) 

The  radially  directed  electric  field,  E^(PN);  i.e. , the  electric  field 
component  normal  to  the  surface  at  P^,  which  is  generated  by  the  source  T 
may  be  obtained  from  Equation  (A-39)  via  the  equation  of  continuity.  Thus, 

7S  * CnxTf^(PN)]  = n • ^(PN)»  (A-40) 


where  [V^.]  is  the  surface  divergence  operator.  The  above  equation  simplifies 
to 
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(A— 41 ) 
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* 
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En(PN)  = " 57<TrHt)+  lz(T2'Rt)  • 

o'—  _ 


which  in  turn  may  be  evaluated  via  Equation  (A-33).  Without  going  through 
the  details  of  the  evaluation  of  Equation  (A-41),  the  final  expression  for 
E^(Pm)  which  contains  terms  only  to  lowest  order  in  ^ is  given  below. 


^(P N>"A>  fcf*1-  fe)3v(0+  (f)  feCu(e)-v(«)] 

- fe(1"  fe-)[“(0-v(«)]ie(i(s). 


(A-42) 


It  is  noted  that  uU)  and  v(c)  are  special  cases  of  the  more  general 
functions  FsU,yj,y2)  and  F^(?,y.|,y2)  which  are  indicated  in  Appendix  V. 

It  can  be  shown  that  if  the  height  of  the  observation  point  P above  the 
surface  is  d2  (i.e.,  jP^P |=d2  of  Figures  1 and  2)  and  d2>0,  then  the 
functions  u(c)  and  v (c)  in  Equations  (A-34),  (A-36b),  (A-39)  and  (A-42) 
may  be  replaced  by  FsU,0,y2)  and  FhU,0,y2),  respectively  with  y-j=0  and 
y2=m"^kd2,  provided  kd2«kpg(PN).  Furthermore,  if  the  source  T of  Equation 
(A-37)  is  also  raised  above  the  surface  from  P^  to  P*  (i.e.,  |P^P'|=d-|  of 
Figures  1 and  2)  when  P^  is  raised  to  P above  the  surface,  then  the  u(0  and 
vU)  in  Equation  (A-39)  and  Equation  (A-42)  must  then  be  replaced  by 
F$(5,y1,y2)  and  Fh(c,y1,y2),  respectively  with  y1=m"1kd1  and  y2=m_1kd2,  for 
kdj«kpg(P|!j)  and  kd2«kpg(PN).  An  approximation  based  on  Taylor  series 
is  given  in  Appendix  V for  F$  foU.y-1.y2)  when  h^  and  h2  are  sufficiently 
small  and  when  #0  (also  5 must  not  be  close  te  zero). 
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APPENDIX  II 

ASYMPTOTIC  ANALYSIS  OF  THE  SURFACE  FIELDS  EXCITED  BY 
INFINITESIMAL  ELECTRIC  OR  MAGNETIC  CURRENT 
SOURCES  ON  A PERFECTLY-CONDUCTING  SPHERE 


The  geometry  of  this  canonical  problem  is  illustrated  in  Figure  3. 
The  case  of  an  infinitesimal  tangential  magnetic  current  source  at  P^  will 
be  treated  first;  a corresponding  treatment  will  be  subsequently  presented 
for  the  case  of  an  infinitesimal  radial  electric  current  source  at  P^. 

The  notation  employed  in  this  Appendix  for  the  fields  E"1,  iff?  due  to 

n x 

the  source  FT  has  the  same  meaning  as  that  in  Section  II,  and  Appendix  I, 
respectively.  Let  the  source  FT  be  x-directed  in  Figure  3;  then  the  vector 
potential  ro  associated  with  the  source  M=xpm5 (r-r1 1 ) in  free  space  (i.e., 
in  the  absence  of  the  sphere)  is  given  by  [14] 


- Pm  0-Jk|r-r' 
F = £ > e_ 

° ^ |- p | 


(A-43a) 


- kPm  “ t?\ 

T0  = Io(2n+l)h^;(kr')3n(kr)Pn(cose);  |r|<lr'|,  (A-43b) 


where  jn,  h^  ' and  Pn  are  the  usual  spherical  Bessel,  Hankel  and  Legendre 
functions.  Let 


r1  = (a+d^ )z  = bz. 


(A-44) 


It  is  convenient  to  introduce  a new  set  of  potentials  trer  and  7rmr  which 
satisfy  the  following  differential  equation  and  the  field  relations  [14]. 
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(A-51) 


where  Er,  EQ,  and  E^  are  the  r,  0 and  <J>  components  of  the  electric  field, 
F which  are  respectively  generated  by  the  potentials  ire  and  irm.  Likewise 
H , H0  and  H are  the  r,  e and  $ components  of  the  magnetic  field,  TT  which 

^ p m q m 

are  respectively  generated  by  it  and  w . Let  n and  V!  denote  the  value 
e m 00 

of  the  potentials  n and  v which  are  associated  with  the  source  FT  in 

free  space.  Let  the  fields  associated  with  these  potentials  iTe  and  tt"1  be 

-fnl  -(n)  00 

denoted  by  E^0'  and  ir0'.  From  [14], 
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H°  = 0 _ 0 o'  0 


(A-52a;A-52b) 


Without  going  through  the  details,  the  following  expressions  are  obtained 
via  (A-43,  A-52a;A-52b). 


r-F^  * rr*F^  = rE^0^  = sind>  -£ 


(A-53) 
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(A-54) 


A useful  procedure  given  by  Fock  in  his  development  of  a modal  expansion 
for  the  fields  of  an  electric  current  moment  over  a dielectric  sphere  [13] 
is  employed  in  order  to  obtain  the  form  in  Equations  (A-53)  and  (A-54). 

The  following  relations  are  helpful  in  arriving  at  Equations  (A-53)  and 
(A-54). 
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Equations  (A-46), 

, (A- 

47),  (A-! 

53)  and  (A-54),  it  can  be 

seen  that 
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the  operator  a* 

is  as  definec 

1 by  Fock 

[13]. 

A = sine  7e  (sine  ^ + ~ J 
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(A-62) 
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p m 

In  order  to  solve  for  TrQ*  in  terms  of  Fq  via  Equations  (A-60)  and  (A-61), 
respectively,  it  is  convenient  to  introduce  the  following  expressions  as 
done  previously  by  Fock  [13]. 


ir*  3 A irj 

T " ■ sin<j>  Te  * T " ■ cos<t>  *Te 


(A-63) ; (A-64) 


Incorporating  Equations  (A-63)  and  (A-64)  into  Equations  (A-60)  and  (A-61), 
with  Fq  as  in  Equation  (A-43b),  it  follows  that 


Ao  ■ - JkY0<TST>  J,  KTHJTT  t'n21(kr,)Jn(kr>pn(cos6) 


(A-65) 
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jn(kr)  Pn(cose)  (A-66) 


where  r*  * b = a+d^. 


Using  the  definitions  [13,15] 
^2)(kr)  = krh<2)(kr)  ; 
one  obtains 


♦n(kr)  = krjn(kr) 


(A-67) ; (A-68) 


Ao  - - JkYo(^r)  wWl  j,  5TWTT  <S2)t«>)*n(kr)Pntcoso)  (A-69) 

Bo  ‘ (l* f)  (kb)tkrT  J,  HT^TT  cnZ>  W>  )♦„(!“•  )pn(<»“)  . (A-70) 

From  the  relationships  in  Equations  (A-63)  and  (A-64),  and  the  results  in 

Equations  (A-69)  and  (A-70),  the  expressions  for  and  ir™  are  completely 

determined.  The  relationship  IL-  P (cose)  = - P^(cose)  may  be  employed  in 
3 A 3B  n n 

evaluating  —■  and  in  Equations  (A-63)  and  (A-64),  respectively. 
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pm 

The  potentials  ir  and  7rQ  are  the  potentials  associated  with  the 
source  FT  in  free  space;  consequently  one  may  Introduce  n®  and  to  denote 
the  potentials  resulting  from  the  presence  of  the  conducting  sphere  such 
that  the  total  potentials  tt®  and  irm  associated  with  FT  in  the  presence  of 
the  sphere  are  obtained  via  the  following  superposition. 


e _ e , e 

IT  ~ 1T_  + IT.  ; 

o s 


m m , m 

IT  = IT  T IT 

0 s 


(A— 71 ); (A-72) 


Let  As  and  B$  be  introduced  as  follows. 


a A.  ^ SB 

T = - sin*  TT  * “F  = “ cos*  se1 


(A-73) ; (A— 74) 


Equations  (A-73)  and  (A-74)  have  been  chosen  to  possess  the  same  form  as 
Equations  (A-63)  and  (A-64),  respectively.  It  follows  from  Equations 
(A-69)  and  (A-70),  that 


kP 


-m.w  / nu  1 v 2n+l 
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(kb)Bn(kr)Pn(cose),  (A-75) 


n _ ,kpnh  k " 2n+l 
Bs  ' ^ fkb)TkrT  nln+TT 


42)’(kb)?(n(kr)Pn(cose) 


(A— 76) 


where  the  coefficients  Bn(kr)  and  An(kr)  must  be  found  by  enforcing  the 

boundary  condition  that  nx[F^+F^]=0  on  the  spherical  surface  at  r=a; 
here,  F^  is  the  electric  field  produced  by  the  potentials  tt®  and  tt™.  The 
total  electric  field  E^F^+F^  is  produced  by  FT  in  the  presence  of  the 
sphere.  Also,  rr*Fn=F^  of  the  earlier  notation.  It  is  easily  verified  that 

n ( ka ) / o \ 

ft(kr)  - ♦ (kr)  - -m — '(kr)  , (A-77) 

" " cj2,(ka)  " 

and 

\ 1>l(ka) 

Bn(kr)  = * (kr)  - -9*p c^'(kr)  . (A-78) 

n n ?(2)  (ka)  n 
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Let  H1”  denote  the  total  magnetic  field  of  fT  in  the  presence  of  the  sphere. 
Then,  from  Equations  (A-63),  (A-64),  (A-65),  CA-66),  (A-71),  (A-72),  (A-73), 
(A-74),  (A-75),  (A-76),  (A-46)  and  (A-49),  it  may  be  seen  that 

jo(2„.1)c(2)(kb)  pj(cose)  (A-79) 

and 


- « kpm  -jkY^cos  <J> 


* /o\i  An(kr)  , 

■]  l (2n+lkJ;  ' (kb) ^PMcose), 


(A-80) 


the  sum  over  n in  Equations  (A-79)  and  (A-80)  starts  from  n=0  instead  of 
n=l  for  convenience  since  P^(cose)=0.  In  order  to  evaluate  Equations 
(A-79)  and  (A-80)  asymptotically  for  large  ka,  it  is  convenient  at  this 
juncture  to  employ  the  Watson  transform  to  these  equations  [5],  As  a first 
step,  one  may  write  Equations  (A-79)  and  (A-80)  more  compactly  for  later 
convenience. 


Figure  A-IV.  Contour  of  integration  in  the  complex  v plane. 

Employing  the  relationship  (-1 )vP^(cose)  = v(v+l )P~^ (-cose),  and  replacing 
v by  -v-1  in  the  integration  over  C*,  one  obtains 


--JE 


dv(iiinr)(2v+1)v(v+1) 


^2)(kb)gfv(kr)/(kr)2(kb) 

\ 


<2)'(kb)K(kr)/(kr)2(kb) 


P”  (-cose). 


(A-84) 


where  e is  a positive  number  however  small.  In  arriving  at  (A-84),  use  is 
made  of  the  following  relationships. 


plJ_l(-cos0)  = P”1 (-cose). 


(A-85) 


• R*  ■£}„<«»  W(W,.  Okr,  . 


(A-86) 


c(2)‘  (kb)X  (kr)  = kb^2)/9(kb)  J , /9(kr)-  Jyi¥-Z~  Hffi/9(kr)  , 


\+v2^0)^/z^r)-  Vi/2^r^"r  * 

(A— 87) 
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2x  l3x  + W * 


X = 


» 


lA-OO 


- 


w = 


F 1 


and  the  circuit  relationships  for  the  cylindrical  Bessel  and  Hankel 

gj 

functions,  Jv+i/2^x^  anc*  ^v+j/2^x^*  resPectively.  Furthermore,  one  may 
employ  the  following  relationships  in  Equation  (A-84);  namely 


8..(kr)/(krr  \ ,2  , , 

- <,—?*  kZ)[-5— ■] 

*2  I dr*  '■* 


&,.(kr) 


A (kr)/(kr) 


kv(v+1)  l (kr) 


V 


(A-89 


P m m 

in  order  to  obtain  integral  representations  for  ir  and  it  , since  En  and 


are  simply  related  to  ire  and  tt111  via  Equations  (A-46)  and  (A-47), 


respectively.  Thus,  from  Equations  (A-46),  (A-47),  (A-81),  (A-84)  and 
(A-89), 


a kpm  -jY  sinij»  . 
e _ / m\/  J o \ 1 

" - — RE — “ 


f dv(5^rK2v-l)£f>(kb)?!v(kr)p;1(-cose) 

--je 

(A-90 


and 


k *-je 

dv(sTnvir)(2v+1)5v2^  (kb)Av(kr)P”1  (-cose). 


m _ / 

IT  = ( 


-»-JE 


(A-91 


irc  is  next  evaluated  for  the  special  case  of  interest;  namely  r'=b=a 


(i.e.,  for  d^=0)  and  r=a  as  follows.  Let  u=v+l/2;  then. 


a kpm  -j  Y si  n 4>  , 
e _ i _n»\#  o \ 1 


:-Je 


_e  _ / m\ / 

* " (^)( — kF 


> d“fe>2“ti-l/2(ka>?!b-l/2(ka)',i-l/2{- 


c 


-r 


i 

r 

ft 

■jt 


i 

f. 

i. 

i- 


In  the  shadow  boundary  transition  region  (l.e. , the  penumbra  region)  and  In 
the  deep  shadow  region,  both  of  which  are  of  interest  in  this  study,  the 
dominant  contribution  to  the  Integral  In  Equation  (A-92)  occurs  when 
lAO(ka).  Thus,  in  this  region,  one  introduces  the  usual  transformation 

Uj 

u=ka+mt  to  approximate  the  cylindrical  functions  ^y_-| /£  and  S-l/2  the 
Fock  type  Airy  functions.  In  particular  [15], 

/T  ka  ^ 

*11-1/2^^  ^ (wi  (t)-W2(t))  ; u=ka+mx  ; m=(7j— ) • (A-93) 

^-l/2(ka)  * - jrp  (W{(t)-W£(t)).  (A- 94) 


(1) 

?y-l/2^ka)  * wi(T)*  (A- 95) 

2 


(1)' 

4-l/2(ka)  ^ 

j/m  2 


(A-96) 


Noting  that  = 2j  e"j(vi"1/2^  £ e-j^'1/2)(27rA)  in  Equation  (A-92), 

one  may  retain  only  the  A=0  term  for  large  ka  as  the  terms  corresponding 
to  contribute  to  the  multiply  encircling  field  paths  around  the  sphere; 
these  terms  (a^O)  contribute  negligibly  for  large  ka.  Also,  Pyli /2 (~cose ) 
may  be  approximated  for  large  ka,  0^tt,  and  iA0(ka)  as  follows  [15], 


P 


-1 

u-1/2 


(-cose) 


7rsin0 


-ju(2ir-e) 


0 


(A-97) 


Let 


ire  = ir®  + ire  (A- 98) 

where  the  term  is  associated  with  the  e-JW0  term  in  Equation  (A-97); 
and  ire  is  associated  with  the  term  e'Jw^2lT"e^  term  in  Equation  (A-97). 
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Simplifying  the  expression  in  Equation  (A-99a)  yields 

*+  v(0  sin$D"  G(ks")  (A-99d) 


where  v(c)  is  defined  in  Equation  (A-29);  it  is  also  discussed  in 
Appendix  V.  Also,  G(ks)  is  defined  in  Equation  (A-27).  Here, 


Clearly  G(ks+)  corresponds  to  the  field  propagation  along  a great  circle 
(geodesic)  ray  path  corresponding  to  s+=ae;  whereas,  G(ks”)  corresponds  to 
propagation  along  the  ray  path  s“=a(27r-e)  along  the  same  great  circle 


(geodesic)  but  in  a direction  opposite  to  s+.  In  this  analysis,  0<e<ir, 


hence,  it  is  noted  that,  D"  = 


sin  (Z-jr-e  J 


2tt— 6 

sin(2ir-e] 


; the  factor  j 


outside  the  radical  is  the  same  as  that  in  front  of  the  term 


in  Equation  (A-97),  and  it  results  from  the  crossing  of  a caustic  of  the 


great  circle  ray  paths  as  they  traverse  the  spherical  surface.  This 
caustic  occurs  at  e=ir.  The  asymptotic  approximation  in  Equation  (A-99)  must 
be  modified  near  6=*;  a simple  modification  of  Equation  (A-99)  for  e at 


and  near  * is  discussed  in  Appendix  III  on  caustic  corrections. 
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Without  going  through  the  details,  it  can  be  similarly  shown  that 


m -pe 
m . rm 

w±'^±-TOi7 — 
o 


tt  £*  r r*  w^t)  — 

6(ks  )D  cosWj—  J dxWgU)  W^-r-yg)-  W2(x-y2)  e_jC  T 


for  r’=b=a  (i.e.,  d^O),  and 
y2=m”^k(r-a)^m”^kd2;  r>a 


The  subscripts  - on  tt+  have  the  same  meaning  as  that  in  tt+;  i.e.. 


m m . m 

IT  = TT+  + IT 


(A-101 ) 


(A-102) 


(A- 103) 


From  Equations  (A-50),  (A-51),  (A-98),  (A-99d),  (A-101),  and  (A-103), 
one  obtains  the  following  surface  magnetic  field  associated  with  tt®  and  tt^ 
(the  field  corresponding  to  it®  and  ir^  is  obtained  by  simply  replacing  s+ 
by  s",  and  by  replacing  £+  by  etc.  in  the  fields  corresponding  to 

_e  ..J  U...  4-U_  J - -J  £ .■ 


H0  = t-H  = Pmcos<j> 


H!J(PN)  * t(t-TT)  + b(b-W). 


subscript  (+)  is  dropped  for  convenience. 

4 

\ 

|i)u(e)*o2  jLvujoG(ks) 

(A-l 04a) 

i 

»[(1-  jj|-)v(C>*.j2D2  ^^|l)G(ks). 

JD 

O 

1 

C 

* 

b(b*W). 

(A-l 04c)  j 

In  obtaining  Equations  (A-l 04a ; A-l 04b ) use  is  made  of  the  following 
relationships. 

a o-Jks"  8 p-Jks*  3 -jks*  i , ± f el 

— + + — - a = 2m3 j-  [-j r];e=<  >.(A-105a) 

36-  ks1  3s”  ks-  ks  ks1  |2tt-0| 


-iy  dVJ{*thi  -4(dVJ{’,)i.  fe  e-Je*,+0*(-iT)e‘J<’T  . 

36-  3 C"  [_3  5 

As  before,  only  terms  to  lowest  order  in  are  retained  in  Equations 
(A-104a;A-104b). 

In  order  to  obtain  E^(PN),  i.e.,  Fr  due  to  FT,  it  is  convenient  to 
rewrite  Equation  (A-46)  as 


(A-l 05b) 


3kV 

It  is  noted  that  a*  is  as  in  Equation  (A-62),  and 


* 2 4 2 2 4 

A _ VT  3 ,S  3 X . D4  3^  3^  + 1 3 

7 - ri? (7 ??> + 7 ;7  ?s?'7 

r=a 


(A-l 06) 


(A-l 07) 


3 ,1 


since  — »•  -*■  -1  (due  to  si n <f>  or  cos<|>  variation  of  ir  * ),  and  (-g-)  is 
3<T  3S  DZ 

neglected  as  it  contributes  to  terms  which  are  higher  order  in  One  may 

rewrite  Equation  (A-106)  via  Equation  (A-l 07 ) as  follows  after  operating 

2 

6 do 

on  ir  with  the  and  — *■  operators. 

9s  3s^ 


-A%e  _ "jpm  I . 2 f , j \ 


-jks  9 . -jks 

+ k2(l-D4)  ^ Dv(5). 

ks  (ks)3 


(A-l 08) 


One  may  conveniently  neglect  the  term  involving  (1-D j in  Equation  (A-108) 

which  contains  a — — T dependence  as  it  is  vanishingly  small  for  (ks)  small 
(ks)3 

(where  D*!)  and  it  is  also  snail  in  comparison  to  the  first  two  terms  for 


(ks)  large.  Thus,  r *f=Er=E1^  becomes 


fri 


Q-jks  I p_ • b 1 = 

^-cr-  • D v(c)(-sin«j>) ; / 
kS  | -si n<^>  P 


(A-l 09) 
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As  before,  the  field  in  Equation  (A-109)  is  associated  with  ?r®  with  the 
subscript  (+)  dropped  for  convenience.  The  term  corresponding  to  ne  is  the 
same  as  in  Equation  (A-109)  except  that  s+  and  s+  are  now  replaced  by  s" 
and  e",  respectively  for  the  it®  case. 


One  may  next  develop  similar  expressions  for  the  electromagnetic 
surface  fields  of  a radial  electric  current  T on  a perfectly-conducting 
sphere.  The  source  7=pe5{r-r')  as  in  Equation  (A-37).  Since  pe=pgn  at  P^, 
the  tangential  magnetic  field  Tl|(PN)  due  to  7 may  be  simply  obtained  from 
r^(PN)  due  to  FT  via  reciprocity  as  indicated  in  Equation  (A-38).  Thus, 

Pg* (n'b[l-  ^-]v(0)D  G(ks).  (A-110) 

One  may  next  compute  the  normal  (radial)  electric  field  ^(PN)  due  to  7 
via  the  equation  of  continuity  (see  Equation  (A-40))  and  Equation  (A-110) 
with  the  understanding  that  b=-<f>  in  Equation  (A-110).  Thus, 

Y[nxHp  *-5 J + r (A-lll ) 


after  neglecting 


3 i 


as  before. 


Then  from  Equation  (A-40), 
1-  fe-)]v(5))D  G(ks). 


(A-112) 


Once  again,  referring  to  Figure  3,  it  is  noted  that  if  d2>0,  then  u(?) 
and  v ( 5)  in  Equations  (A-104),  (A-109),  (A-110)  and  (A-112)  must  be  replaced 
by  Fs(5,0,y2)  and  F^(c,0,y2)  of  Equations  (A-139)  and  (A-138)  for  kd2«ka. 

On  the  other  hand,  if  d^>0  and  d2>0  in  Equation  (A-110)  and  Equation  (A-112); 
then  one  must  replace  u(s)  by  Fs(5,y-|,y2)  and  v(5)  by  Fh(5,y-|,y2)  in  these 
equations  when  kd^  2«ka.  The  functions  Fg  ^U,y-|,y2)  are  discussed  in 
Appendix  V along  with  an  approximation  for  these  functions  based  on  a 
Taylor  expansion  which  is  valid  for  kd-j  2 sufficiently  small  and  cfO. 

Finally,  it  is  once  again  noted  that  the  results  in  Equations  (A-104), 

(A-109),  (A-110)  and  (A-112)  are  not  valid  at  and  near  e=ir.  For  e-*ir, 

these  results  must  be  modified  as  shown  in  Appendix  III  on  caustic  corrections. 
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APPENDIX  III 

CAUSTIC  CORRECTIONS  FOR  THE  SURFACE  FIELDS  ON 
THE  SPHERE 


As  pointed  out  at  the  end  of  Appendix  II,  the  expressions  for  the 

surface  field  on  a sphere  due  to  a source  at  P^  on  the  same  sphere 

(see  Figure  3)  as  given  in  Equations  (A-104),  (A-109),  (A-110)  and  (A-112) 

are  not  valid  at  e=ir  because  D+*  as  e-^ir.  The  point  P^  on  the  surface 

corresponding  to  e=ir  is  a caustic  of  the  surface  rays  since  at  e=ir,  every 

great  circle  path  is  now  a geodesic  between  the  diametrically  opposite 

source  and  field  points,  and  P^,  respectively.  Such  a continuum  of 

geodesic  surface  ray  paths  intersect  at  P^,  thereby  producing  a point 

caustic  at  for  e=ir.  In  order  to  obtain  an  asymptotic  approximation 

for  the  surface  field  which  is  valid  at  and  near  the  caustic  at  e=ir,  one 

begins  by  noting  that  the  approximation  in  Equation  (A-97)  is  responsible 

for  the  singular  behavior  of  D at  6=ir.  Thus,  to  remedy  this  singular  behavior 

of  D at  e=ir,  one  may  employ  an  asymptotic  approximation  for  the  Legendre 

function  which  is  different  from  that  in  Equation  (A-97)  when  e-hr;  such 

an  approximation  is  indicated  in  [15],  Alternatively,  one  may  employ  an 

approximation  in  Equations  (A-104),  (A-109),  (A-110)  and  (A-112)  for  e+ir 

which  is  based  on  the  following  physical  considerations.  Firstly, 

whenever  e^ir,  there  are  two  great  circle  geodesic  surface  ray  paths  between 

0 m 

the  points  P^  and  P^;  these  paths  are  associated  with  the  potentials  tt+* 
and  respectively  as  indicated  earlier  in  Appendix  II.  Furthermore, 

the  total  field  is  constructed  from  the  sum  of  these  tt®*111  potentials  as 
indicated  in  Appendix  II,  i.e.,  from  7re,m  with 


Tre,m  = Tr;*m  + ir®*m  . (A-113) 

0 m 

The  electromagnetic  surface  fields  resulting  from  it  * are  given  in 
Equation  (A-104)  or  (11)  and  in  Equation  (A-109)  or  (15)  for  the  case  of  the 
magnetic  current  source  FT  on  a sphere;  and  they  are  given  in  Equation  (A-110) 
or  (22)  and  Equation  (A-112)  or  (24)  for  the  electric  current  source  J on  a 
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sphere.  From  (11)  it  is  noted  that  the  total  El?  field,  (i.e. , including 
the  ir_’m  contribution)  contains  the  following  combination  of  D G(ks)  type 
terms. 

Tcl  = D+  G(ks+)  + 0“  G(ks"),  (A-114) 
and 


Tc3  = [D+]3G(ks+)  + [D"]3G(ks"). 


± + 


Simplifying  T ^ and  Tcg  via  G(ks  )D 


(A-115) 


(A-116) 


(A-117) 


with  the  understanding  that: 

0+=e  ; s+=ae  ; q~=Zv-q  ; s"=a(2ir-e).  (A- 1 1 8a ; A- 1 1 8b ; 

A- 1 1 8c ; A- 1 1 8d ) 

Since  0<e<ir  is  understood  in  the  present  development,  it  is  clear  that 
s1ne~<0.  Furthermore, 


e" 

2*-e  _ - , [ 

'sine”  ' 

-sine  J \| 

2ir— 8 
sine  ’ 


O<0<TT 


(A-119) 


hear  the  caustic,  e-*ir ; hence,  Equations  (A-116)  and  (A-117)  may  be 

epproxlawted  as 


kZyo  H-  r-.1kaw^ka<’,-e)  , e-jka(-9)\l  , 

**rL  v*s~  * 3 “S5- JJ  17 1 


T = J- 
c3  sine 


Tc2  ’ 


(A-120) 


(A-121 ) 


where 


c2  * 57T 


x r.  c,ik»pa<*-9>  . e-Jka t"'9 )Ti , . 
*j  If 6 V /JTe  J ^5-  /jks+  • 


e+TT  . (A- 122) 


In  the  above  approximations,  sine  = sin(Tr-e)^  tt-g  has  been  used  for  e-*ir; 

also  ^—7%  — r is  used  for  e-nr.  One  may  now  rewrite  Equation  (A-120)  and 
ks~  ks 

(A-122)  in  terms  of  trigonometric  functions  as: 

Tcl  ^ 77,113/2  e_j2Trni  “?  (A‘ 


(A-123) 


Tc2  * TT  7rm3/2  e"j2lirT1  sin(ka(w-e)-  J)l  . 


(A-124) 


kp  1/3 

m = (-^r1)  , and  p =a  for  the  sphere. 


(A-125) 


It  is  noted  that  sin(ka(ir-e)-  ^ ) = cos(ka(it-e)-  \ - *r)  in  Equation  (A-124). 
Furthermore,  one  also  notes  that  for  ka (ir-e )>>  1 the  trigonometric  functions 
in  Equations  (A-123)  and  (A-124)  are  asymptotic  approximations  of  cylindrical 
Bessel  functions  as  indicated  below. 
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I 

i 


J0(kaU-e))  * 


2 

rka  (tt-0  ) 


cos(ka(ir-e)-  J)! 


(A-126) 


and 


J] (ka (tt-0 ) ) _e ) cos(kaU-e)-  J- - £-)> 


(A-127) 


For  small  kaCn-e),  one  may  therefore  replace  the  trigonometric  functions 
in  Equations  (A-123)  and  (A-124)  by  the  Bessel  functions  in  Equations 
(A-126)  and  (A-127),  respectively,  thereby  making  and  Tc2  (and  hence 
Tc3)  va^1c*  ^or  This  Procedure  is  similar  to  that  employed  earlier  by 

Keller  [3]  in  his  study  of  the  diffraction  by  a circular  aperture. 

Thus,  for  e^ir,  one  may  employ  the  following  expressions  for  T J, 

Tc2  and  Tc3  instead  of  those  in  Equations  (A-120),  (A-122)  and  (A-121), 
respectively. 


From  Equations  (15)  and  (22)  it  is  noted  that  the  total  ^ and  the 
fields  (i.e.,  including  the  ir®,ni  contribution)  respectively  contain  only 
the  type  term;  whereas,  the  total  ^ field  corresponding  to  Equation 
(24)  contains  only  the  Tc-|  type  term.  Hence,  Equations  (15)  and  (22)  may  be 
modified  for  e-*ir  by  employing  of  Equation  (A-129)  in  place  of 
Equation  (A-122);  likewise  Equation  (24)  may  be  modified  for  e-*-ir  by  employing 
Tcl  of  Equation  (A-128)  in  place  of  Equation  (A-120). 
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APPENDIX  IV 

NUMERICAL  TECHNIQUES  FOR  EVALUATING  THE  MULTIPLE 
SURFACE  INTEGRALS  ASSOCIATED  WITH  THE  MUTUAL 
ADMITTANCE  BETWEEN  TWO  IDENTICAL  SLOTS 
ON  A CIRCULAR  CYLINDER 


l 

I 


Consider  two  identical  circumferential  slots  on  the  surface  of  a 
circular  cylinder.  Under  the  "dominant  mode"  approximation  of  the  aperture 
fields  in  the  slots,  the  mutual  admittance  between  the  two  slots  can  be 
expressed  as  follows: 


,c  -2a2 
12  “ id 


J d<h 

<t>,=-t,/2a 


dz. 


z,=-d/2 


4>  +£/ 2a 
Or 


di|>. 


zo^2 


dz. 


*2=V*/2a  z2=Zo"d/2 


ira<j>. 

cos[— cos[|a(*2-<f>0 


)]  gj(s.e) 


(A-131) 


In  the  above  expression,  a is  the  radius  of  the  circular  cylinder,  d and  a 
are  the  slot  dimensions;  and  <|>0  and  zQ  are  the  angular  and  axial  separation 
of  the  centers  of  the  two  slots,  respectively.  These  parameters  are 
illustrated  in  Figure  A-5.  The  Green's  function  g£(s,e)  represents  the 
surface  field  at  (<f>2*z2 ^ due  to  a unit-strength,  ^-directed 

magnetic  dipole  at  UpZ^).  The  variable  "s"  is  the  arc  length  along  the 
geodesic  path  joining  the  source  location  (^pZ-j ) and  the  field  point 
location  (<|>2,z2);  and  the  variable  "0"  is  the  pitch  angle  associated  with 
the  geodesic  helix  as  shown  in  Figure  A-5.  It  has  been  shown  in  Section  II 
that  the  surface  field  is  indeed  only  a function  of  the  arc  length 

"s"  and  the  pitch  angle  "e".  It  is  noted  that  both  "s"  and  "e"  depend 
upon  the  relative  distance  of  jz2-z-j[  anc*  1 4>2_<t’i  l * The  feature  that 

depends  upon  only  the  relative  axial  and  angular  separation  between 
the  source  and  field  locations  suggests  an  efficient  way  to  perform  the 
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X 


Figure  A-V.  Slot  dimensions  and  the  cylinder  geometry. 

double  surface  integration  required  for  the  evaluation  of  the  mutual 
coupling  between  the  slots.  In  fact,  the  double  surface  integration  can  be 
reduced  to  a double  line  integration  via  a suitable  coordinate  transformation. 
Equation  (A-131)  can  be  rewritten  as  follows: 


2 4 


^ l l' 


(A-l 32a ) 


where 


I-j  = cos(~  $Q)  j j ||  cos[^-U2+^1)]gJ(s,e)dz1d22d^1d^2  (A-132b) 

*1*2Z1Z2 


(A-1 32c) 


n 


l2  s cos(r-  *o)  J { J J cos[|^(^2-<))1)]gJ(s,e)d21dz2d<t>1d^ 

V2Z1Z2 


I3  = sin(~  <J>0)  Jfff  sint—C^^lgJCs.eJd z1dz2d#1d*2  (A-132d) 

<t*1<t>2zlz2 


u = sin(^^0) 


////  sin[p<^2-<(>1  )]g^(s,e  )dz1dz2d<j>1d<j>2  (A-132e) 


Kr2zlz2 


and 


s =>|  a (<|>2-<t>1 ) + (z^) 


e = tan 


-1 


Z2-Zl 


(A-1 33a) 
(A-1 33b) 


It  is  observed  that  all  four  integrals  in  Equations  (A-1 32b; A-1 32c; A-1 32d; 
A-132e)  are  similar  in  form.  This  appendix  is  devoted  to  the  discussion 
of  the  numerical  techniques  used  to  perform  only  the  double  surface 
integration  appearing  in  Equation  (A-132b).  Similar  analysis  can  be 
carried  out  for  the  rest  of  the  integrals  in  Equations  (A-132c;A-132d;A-132e). 


Let  us  introduce  the  following  coordinate  transformations. 


— - ) 

(A-1 34a) 

” ($2+<h ) 

(A-1 34b) 

X (Z2’Z1 1 

(A-1 34c) 
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J 


(A-1 34d) 


w = — (z?+z, ) 

SZ  c 1 

It  can  be  shown  that  the  integral  I-j  given  in  Equation  (A-l 32b),  via  the 
transformation  defined  in  Equations  (A-l 34a ; A-l 34b; A-l 34c; A-l 34d),  can  be 
expressed  as  follows: 


Vo  U9(t)  t~  Wp(t) 


Ij  = cos 


v2  U2VW  2 2V  ; 

(f^0)  J J J f cos[pV?u]gJ(v,t)dtdw  dudv. 

v]  u] (t)  t1  Wl(t) 


(A-l 35) 


The  regions  of  integration  in  the  (t,w)  plane  and  (u,v)  plane  correspond 
to  the  shaded  areas  illustrated  in  Figure  A-6.  The  u-integration  and 
w-integration  in  Equation  (A-l 35)  can  be  readily  integrated  in  closed 
form  and  leads  to  the  following  expression. 


v?=  — U +*/a) 

r sz  0 

\ dv 

v ^v*/a) 


‘2-  -^(vd> 


r /2 

J dt  F(t)[7(u2)-F(u1')]gJ(v,t) 


t,=  — -(z  -d) 


(A-l 36a) 


where 


F(t)  = 


2t  - /7(z  -d)  ; t < — 

o 


JZ  (z  +d)-2t  ; t > — 

° “ /Z 


(A-l 36b) 


IT"^= 

? 

* 

t2( 

Ji 

w-t 

V/  V 

w,(t) 


W-t  * yfz  (i.) 


= y/Z  (2e  + 


/I 


\ W + t = y/Z  (20-|-) 


Figure  A-VIa.  Region  of  integration  in  the  (t,w)  plane. 


v*t 


Vz  <4> 


u ,(vr 


-u-v  = 72  (^) 


^2<V) 


/I 


— -u+v  8 -JZ  ( - 


Figure  A-VIb,  Region  of  integration  in  the  (u,v)  plane. 


i»|o. 


F(u)  = 


sln(^f-u) 


(A-136c) 


U1  = 


v - rt  fc/2a 


^(4»0-Jt/2a)-v 


v > 


rt 


*0 

V < — 


(A-136d) 


( */2’(  4>0+*/2a ) -v 


U2  = < 


v + /7U/2a) 


v > 


v < 


rt 


From  Equations  (A-132b)  and  (A-135),  one  concludes  that 


(A-136e) 


1/ 2a  d/2 


dz. 


<0o+t/2a 


d<f>5 


'=-a/2a  z^— d/2  4»2^4>0"*’/2a 


z +d/2 


dz2  cos[  j-(<f>2+<(.1 ) ]g^  (s  ,e ) 

z^Zo-d/2 


1 


■— (a  +£/a)  ~(z  +d) 


rtf0  ' /?° 

J dv  j 

V ~(^0-^/a)  t=  — < 
rt  0 


dt  F(t)[F(u2)-F(u1)]gJ(v,t)  (A-137) 


The  above  equation  indicates  that  the  double  surface  integration  in  Equation 
(A-l 32b ) can  be  reduced  to  a double  line  integration  via  a set  of 
suitable  coordinate  transformations. 


It  should  be  noted  that  the  discussion  in  this  appendix  is  also 
valid  for  the  case  of  two  axial  slots  on  the  circular  cylinder. 
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APPENDIX  V 

ON  THE  SURFACE  TYPE  FOCK  FUNCTIONS 


In  this  appendix  we  define  and  list  some  useful  formulas  of  the  Fock 
functions:  FhU,ylty2).  FsU,ylty2),  vU),  uU),  v7U).  v'U),  u'U),  and 
the  Airy  functions  w-j(t)  and  w2(t).  For  a complex  x and  real  e,  we  define 
the  above  functions  as  follows. 


FhU,yi,y2)  = 5j 


L L_  ej*  51/2  • 


/3tT 


r .Wi(t)  ■ _ 

• J dx  WgCx-y^CW^x-yg)  - W2(x-y2)]e"J  T 

r7  2 l 2 1 

yl  > y2 
2 1 


Fs(C’yly2^ 


rrr  „ o 


(A-138) 


y 1 /_  w W,(x) 

yf> 


■j^T 


yl  ” y2 
2 1 


where  y7  and  y2  are  defined  in  Sections  II  and  III.  For  y7*y2=0: 
F„(e,0,0)  ■ v(5)  - J-  e*  51/2  / d,  e-^T 


(A-139) 


(A-140) 


39 


Wo(t) 


FsU,0,0) 


^ 53/2  ^ J dT  4^ 


(A-141 ) 


Vl(5) 


h I 


W2(t) 

wJTTT 


2VW  -iFr 

dx  HTT-r  T e j5t 


(A-142) 


1 f z3 

w^x)  = — I dz  exp{xz-  j-) 

* 

t r z3 

w2(x)  = — J dz  exp(xz-  3-) 


CA-143) 


(A-144) 


where  integration  contour  r^(r2)  goes  from  ® to  0 along  the  line 
arg(z)  = - j^-(+2tt/3)  and  from  0 to  « along  the  real  axis. 

For  &0,  and  y1  ,y2  small,  one  may  approximate  Fg  and  of  Equations 
(A-139)  and  (A-138)  by  a Taylor  series  expansion  as  follows. 


Fh(5.ylty2)  % v(c)  - j r1  v^ej^yf];  C/0 


Fs(£,y1,y2)  * uU)  + 1 [u'U)-  | 5_1uU)3  [yf+yg];  5^0 


(A-145) 


(A-146) 


For  5 large  and  positive,  one  may  employ  a rapidly  converging  residue  series 
representation  for  the  Fock  functions  as  follows  [6,11,12] 


v (c ) = e"J2rv751/2  l (x')"1  e"j5T" 

n-1  n 


(A-147) 
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CO 


(A-148) 


u(c)  - kVZ  l e’JCTn 

n=l 


Jt  0/0  “ “j^Tn 

v,  (C)  « e * 2A  s3/2  l e 
1 n=l 


(A-149) 


j71  ^ j | 

v*  (c)  = \e  4 ^rr1/2  I (1-J2cti;)(t|;)-1  e"  n (A-150) 

n=l 

u'(5)  = e^3v751/2  jf  n)e‘J5Tn  (A-151) 

n=l  J n 

where  xn  and  x^  are  zeros  of  w2(x)  and  w2(x),  respectively,  and  they  are 
tabulated  in  Table  A-I. 


TABLE  A-I 

Zeros  of  W2(x)  and  W2(x) 


x = |xl  e"J>/3  and  x'  = |x'|  e'jir/3 
n ' n ' n n 

- r-  | rrn 


Tn 

lTJ 

1 

2.33811 

1.01879 

2 

4.08795 

3.24819 

3 

5.52056 

4.82010 

4 

6.78661 

6.16331 

5 

7.94413 

7.37218 

6 

9.02265 

8.48849 

7 

10.0402 

9.53545 

8 

11.0085 

10.5277 

9 

11.9300 

11.4751 

10 

12.8288 

12.3848 

On  the  other  hand,  if  £ is  small  and  positive,  one  may  employ  a small 
argument  asymptotic  expansion  for  the  Fock  functions  as  follows  [6,11,12]. 


v(cH-  2T  ej*«3/2+  Q 

i A 7/7  9/2 

r5  + ?TTe  ? 

(A-152) 

uUKi-  ^ e\3/2+  ^ 

L3a  5/S"  "A  9/ 2 

K + p-  e 5 -••• 

(A-153) 

^7  eJ^53/2-  \ 

(A- 154) 

V'(0  ^e'J^51/2+  l 

(A- 155) 

u'uy  l^  e J_2rc1/2+  ] 

jj.2 , 45/7  “jTr7/2 

n + we  ? — • 

(A-156) 

For  ZH.q,  the  residue  series  representation  with  the  first  ten  terms  in  the 
summation  may  be  used.  For  ?<£Q,  the  small  argument  asymptotic  expression 
with  the  first  three  terms  may  be  used.  In  the  present  study  is  set 
to  be  0.6  as  shown  in  [11,12], 
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